Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  809 1859 (N44E53 L=321)   871     69        1 ( 782/1832- 842/1888)  -8.1 E60-E46 N41-N46
      2  545 1761 (N29E64 L=310)  1547    158        4 ( 483/1711- 603/1802)  83.7 E71-E59 N24-N33
      3  415 1606 (N16E62 L=312)  1361    146        3 ( 352/1549- 458/1644) -11.3 E75-E55 N12-N17
      4  916 1570 (N26E22 L=352)   625    143        1 ( 894/1500- 950/1629) -73.3 E24-E19 N20-N30
      5 1272 1232 (N14W09 L= 23)   249    102        0 (1223/1203-1316/1261)   0.0 W05-W12 N13-N15
      6  663  935 (S18E19 L=355)   365    113        0 ( 631/ 912- 705/ 970)   0.0 E22-E16 S19-S16
      7  801  705 (S28E05 L=  9)   762    292        0 ( 659/ 673- 926/ 741)  33.7 E14-W01 S33-S22
      8 1455  411 (S27W50 L= 64)   683     39        0 (1429/ 382-1487/ 437)   0.0 W46-W53 S28-S24
      9 1292  327 (S37W45 L= 59)  2897    370        5 (1130/ 267-1428/ 414)  50.9 W36-W52 S46-S28
  

    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