Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  794 1874 (N55E55 L=212)   665     80        1 ( 762/1845- 831/1900)  -9.9 E69-E45 N52-N56
      2 1047 1758 (N50E16 L=250)   442     60        0 (1018/1735-1077/1778)  68.2 E18-E13 N47-N52
      3 1611 1686 (N58W53 L=320)   739     62        2 (1582/1661-1645/1703)   9.5 W45-W69 N54-N61
      4 1190 1661 (N44E00 L=266)   525     83        0 (1146/1637-1226/1683)   9.5 E03-W03 N42-N46
      5 1672 1377 (N33W43 L=309)   315     30        0 (1649/1364-1695/1393)   0.0 W40-W45 N31-N33
      6 1425 1284 (N22W22 L=288)   477    167        0 (1366/1228-1489/1333) -29.1 W17-W26 N19-N24
      7  175  967 (S17E65 L=201)   634     53        2 ( 155/ 935- 202/ 992)  53.1 E68-E60 S19-S15
      8 1702  681 (S10W54 L=321)   327    101        0 (1679/ 638-1726/ 723)  90.0 W52-W55 S12-S06
      9 1293  478 (S30W30 L=296)   569    130        0 (1225/ 459-1361/ 498)  24.0 W25-W34 S31-S27
  

    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