Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1054 1759 (N52W32 L=290)   613     83        2 (1018/1723-1077/1794)  11.3 W27-W37 N48-N54
      2 1456 1724 (N35W60 L=318)   500     77        1 (1412/1711-1499/1736) -20.6 W55-W64 N33-N36
      3  563 1365 (N38E24 L=234)  1062    221        2 ( 477/1302- 657/1399) -31.0 E31-E16 N32-N42
      4 1063 1347 (N25W11 L=269)   979    336        0 ( 934/1324-1235/1376) -14.7 W02-W21 N20-N27
      5 1776 1146 (S06W53 L=310)   976    120        4 (1742/1088-1803/1209)  45.0 W48-W57 S09-S02
      6  588 1087 (N20E25 L=233)   846    240        0 ( 495/1005- 676/1141) -50.2 E31-E20 N13-N25
      7  208  968 (N20E61 L=197)  1033    113        1 ( 171/ 916- 243/1013) -36.9 E66-E57 N16-N22
      8 1690  642 (S35W37 L=295)  6435    880       12 (1476/ 392-1865/1001)  34.9 W16-W59 S49-S17
      9 1259  877 (S07W10 L=268)   235    102        0 (1212/ 866-1310/ 888) -31.0 W07-W12 S08-S05
     10 1275  662 (S20W06 L=264)   650    278        1 (1189/ 562-1313/ 776)  47.7 E01-W09 S24-S13
     11 1255  262 (S50E08 L=250)  1952    367        5 (1096/ 187-1384/ 356)   9.2 E25-W06 S53-S45
     12  961  248 (S41E30 L=228)   535     65        1 ( 923/ 228- 995/ 275)   0.0 E34-E26 S43-S39
  

    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