Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1169 1737 (N44W06 L= 63)   310     52        0 (1144/1716-1204/1755)  21.8 W03-W08 N41-N45
      2 1284 1679 (N40W15 L= 72)   265     72        0 (1268/1649-1301/1705)  68.2 W13-W16 N37-N42
      3 1367 1666 (N40W22 L= 79)   509     66        1 (1344/1632-1396/1690)  51.3 W19-W24 N37-N42
      4  926 1620 (N32E11 L= 46)   747    157        3 ( 869/1557- 980/1671) -41.6 E15-E06 N28-N36
      5 1467 1623 (N38W29 L= 86)   611    113        0 (1441/1576-1499/1656)  71.6 W27-W32 N33-N40
      6 1644 1519 (N32W42 L= 99)  2743    495        8 (1509/1435-1804/1594)  -9.1 W31-W56 N27-N36
      7  633 1486 (N22E29 L= 28)   244     81        0 ( 616/1443- 648/1522) -59.0 E30-E27 N18-N23
      8  203 1301 (N09E62 L=355)  1269    366        4 ( 147/1226- 254/1378) -32.0 E73-E55 N04-N14
      9  803 1305 (N11E15 L= 42)   239     54        1 ( 770/1284- 830/1325)   0.0 E17-E13 N09-N12
     10  454 1193 (N03E37 L= 20)   243     73        0 ( 422/1168- 488/1217)  36.9 E39-E35 N01-N04
     11  856  696 (S26E09 L= 48)   253     72        0 ( 829/ 660- 884/ 734)  53.1 E10-E07 S27-S23
     12  412  573 (S35E44 L= 13)   845    179        2 ( 351/ 543- 489/ 614) -15.3 E48-E37 S37-S32
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y