Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  592 1661 (N19E48 L=224)   423     52        0 ( 562/1638- 616/1687)  63.4 E49-E46 N17-N21
      2 1241 1623 (N35E05 L=267)   243     42        1 (1228/1596-1258/1647) -56.3 E07-E04 N33-N37
      3  377 1415 (S00E53 L=219)   372     59        0 ( 358/1382- 392/1448) -56.3 E54-E50 S02-N01
      4  855 1357 (N07E19 L=254)   488    153        0 ( 797/1307- 924/1415) -16.7 E23-E13 N06-N09
      5 1320 1206 (N12W11 L=284)  1108    311        1 (1283/1037-1345/1349) -62.1 W05-W16 N02-N19
      6  136 1202 (S16E70 L=202)   676     97        2 ( 118/1174- 154/1230)  63.4 E72-E67 S18-S14
      7  212 1183 (S17E60 L=212)   749    142        2 ( 174/1135- 248/1225)  71.6 E62-E56 S20-S14
      8 1524  936 (N03W31 L=303)   363    117        0 (1491/ 884-1541/ 997) -80.5 W28-W32 N00-N06
      9  356  915 (S30E41 L=232)   973    284        0 ( 258/ 883- 455/ 956)  20.6 E49-E33 S32-S24
     10 1013  872 (S15W04 L=276)  1179    340        1 ( 835/ 830-1180/ 916)  16.7 E07-W13 S18-S11
     11  218  679 (S48E55 L=217)   568     98        0 ( 195/ 664- 252/ 705)  -9.5 E58-E51 S49-S46
  

    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