Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1605 (N37E13 L=355)   339    118        1 ( 928/1563-1022/1646) -29.1 E17-E08 N34-N39
      2 1763 1441 (N37W59 L= 67)   537     58        0 (1735/1422-1797/1462)  -6.3 W54-W64 N36-N38
      3  541 1368 (N15E37 L=331)   579    232        1 ( 503/1246- 566/1457) -55.3 E42-E32 N07-N20
      4  262 1370 (N11E63 L=305)   776     88        1 ( 240/1334- 292/1416) -71.6 E67-E59 N08-N14
      5 1696 1333 (N29W46 L= 54)   537     98        2 (1656/1288-1731/1374) -33.7 W42-W49 N25-N31
      6  790 1199 (N08E17 L=351)   885    276        1 ( 739/1079- 837/1314)  82.4 E18-E14 N00-N15
      7 1018 1171 (N09E02 L=  6)   234    108        0 ( 995/1118-1040/1218) -51.3 E04-W00 N06-N11
      8  474 1105 (S02E36 L=332)   352     70        0 ( 446/1080- 508/1137) -21.8 E38-E33 S03-N00
      9 1008  653 (S23W04 L= 12)  2066    665        1 ( 754/ 555-1287/ 773)  34.1 E12-W19 S32-S11
     10  351  599 (S36E49 L=319)   855    197        0 ( 310/ 524- 400/ 685) -74.7 E52-E45 S41-S30
     11  521  621 (S32E31 L=337)   414     86        1 ( 477/ 608- 564/ 635)   0.0 E34-E27 S32-S30
  

    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