Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1325 1617 (N37W03 L= 24)   307     18        0 (1303/1598-1350/1636)   0.0 W00-W04 N35-N38
      2  902 1421 (N12E17 L=  5)   294     96        1 ( 858/1391- 940/1445)  59.0 E18-E15 N09-N14
      3  323 1117 (S17E47 L=334)  1308    278        1 ( 268/ 991- 358/1258)  86.6 E50-E45 S25-S08
      4  115 1131 (S18E70 L=311)   853    235        1 (  82/1104- 153/1153)  20.6 E74-E65 S19-S16
      5  601 1131 (S11E27 L=354)   361     73        2 ( 558/1115- 639/1153)  36.9 E30-E24 S12-S08
      6  749 1081 (S11E17 L=  5)   229     97        0 ( 707/1037- 785/1104) -18.4 E19-E13 S12-S09
      7 1899 1010 (N18W61 L= 82)  1254    212        4 (1870/ 963-1935/1058) -51.3 W57-W65 N15-N20
      8  838  980 (S14E09 L= 12)   187     77        0 ( 795/ 960- 880/1005)  45.0 E11-E07 S16-S11
      9 1075  831 (S17W08 L= 29)   375    156        1 (1045/ 747-1109/ 908) -48.8 W04-W12 S20-S12
     10 1669  781 (S03W46 L= 68)  1297    263        1 (1592/ 734-1757/ 861)   0.0 W37-W55 S04-N00
     11  756  720 (S31E09 L= 13)   900    297        1 ( 613/ 682- 901/ 753)  11.9 E18-W01 S34-S28
     12 1075  713 (S23W11 L= 33)   246     58        0 (1039/ 688-1110/ 739)  53.1 W09-W12 S25-S20
     13 1271  698 (S19W23 L= 45)   205     73        0 (1234/ 672-1297/ 728)  78.7 W21-W24 S21-S16
     14 1512  453 (S25W50 L= 71)  2587    428        4 (1343/ 383-1712/ 557)  42.1 W40-W63 S33-S14
     15 1213  531 (S30W26 L= 48)   250     63        0 (1182/ 515-1245/ 548)  26.6 W24-W28 S30-S28
     16 1620  391 (S24W63 L= 84)  1639    202        7 (1584/ 336-1653/ 454) -54.5 W59-W66 S27-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