Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  975 1665 (N43W17 L=206)   554    128        1 ( 939/1620-1036/1701)  16.7 W12-W22 N40-N44
      2  430 1304 (N34E34 L=155)  1135    190        4 ( 372/1237- 508/1366)  22.6 E39-E26 N30-N36
      3  983 1286 (N19W04 L=193)   267     86        0 ( 969/1245- 998/1331)  45.0 W02-W06 N16-N21
      4 1877 1034 (S18W60 L=249)  3521    492        5 (1810/ 901-1950/1238) -70.8 W55-W73 S27-S04
      5  796  852 (S01E17 L=172)   297    131        0 ( 780/ 795- 814/ 910)  63.4 E18-E14 S04-N02
      6 1502  766 (S23W23 L=212)   287     77        0 (1484/ 724-1516/ 804)  78.7 W21-W23 S25-S20
      7  192  758 (N06E66 L=123)   532     61        0 ( 170/ 736- 222/ 788) -26.6 E69-E62 N04-N08
      8 1682  688 (S33W36 L=225)  1051    178        0 (1637/ 615-1723/ 769)  48.8 W32-W40 S37-S28
      9  903  641 (S15E16 L=173)   230     82        0 ( 876/ 614- 920/ 680)  45.0 E18-E14 S16-S12
     10 1294  575 (S29W05 L=194)   288     88        2 (1253/ 544-1328/ 606)   9.5 W01-W07 S30-S27
     11  649  346 (S28E46 L=143)   437    120        0 ( 632/ 290- 668/ 400)  37.9 E50-E40 S31-S24
  

    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