Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1553 (N43E03 L=156)   359     90        0 ( 813/1529- 905/1575) -33.7 E07-W01 N40-N45
      2 1731 1480 (N20W64 L=224)  3711    471       10 (1631/1411-1830/1535) -20.0 W58-W73 N13-N25
      3 1333 1316 (N19W25 L=184)   772    205        3 (1263/1248-1415/1354)   0.0 W21-W28 N13-N22
      4 1594 1290 (N12W42 L=202)   321    107        2 (1545/1270-1630/1318)   0.0 W40-W45 N10-N15
      5  358 1271 (N30E47 L=113)   369     90        0 ( 323/1257- 393/1285)  -8.1 E51-E43 N29-N31
      6 1727 1175 (N03W50 L=210)   328    120        1 (1687/1151-1752/1217)  15.9 W48-W52 N01-N06
      7  760 1143 (N17E15 L=145)   596     99        1 ( 731/1083- 787/1198) -90.0 E16-E14 N13-N21
      8 1240 1166 (N11W15 L=175)   309     78        2 (1202/1134-1280/1196) -36.9 W14-W17 N09-N14
      9 1837  956 (S13W59 L=218)  2170    511        9 (1793/ 773-1898/1125)  55.6 W56-W75 S27-S02
     10 1234 1038 (N04W13 L=172)   246     86        0 (1189/1022-1269/1061)   0.0 W10-W15 N02-N06
     11 1265  912 (S04W13 L=172)   504    115        1 (1202/ 889-1316/ 944)   0.0 W09-W16 S06-S01
     12 1348  766 (S14W16 L=176)  1107    227        4 (1310/ 670-1386/ 881)  69.0 W15-W18 S21-S07
     13 1683  773 (S21W43 L=202)   775    147        1 (1645/ 743-1719/ 803)   8.1 W39-W47 S24-S18
     14  923  535 (S23E15 L=144)   554    179        0 ( 905/ 456- 942/ 614)  70.0 E14-E16 S28-S17
     15 1203  502 (S30W04 L=163)  1028    185        0 (1111/ 449-1266/ 573)  24.8 E02-W08 S35-S23
     16  463  450 (S23E55 L=105)   956    184        4 ( 414/ 416- 506/ 494)  12.1 E57-E53 S26-S20
  

    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