Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1578 1697 (N35W58 L=131)   730     67        1 (1541/1670-1609/1727)   8.7 W52-W65 N33-N37
      2  748 1532 (N34E13 L= 60)   261     39        0 ( 720/1511- 774/1549) -36.9 E15-E11 N32-N35
      3  525 1394 (N28E30 L= 43)  1270    383        2 ( 331/1344- 678/1429) -13.1 E48-E18 N22-N31
      4 1718 1297 (N08W50 L=123)   308     53        1 (1705/1263-1737/1328)  76.0 W48-W51 N05-N09
      5  138 1219 (N21E69 L=  4)  1404    106        2 ( 102/1169- 178/1262) -18.4 E82-E61 N17-N24
      6  446 1169 (N14E36 L= 37)   689    209        0 ( 391/1106- 525/1251)  35.0 E40-E29 N10-N18
      7  463 1012 (N05E35 L= 38)   517    181        1 ( 409/ 966- 532/1063) -32.0 E39-E30 N01-N08
      8  561  946 (S00E29 L= 44)   433    171        0 ( 474/ 924- 623/ 989)   5.2 E35-E24 S01-N01
      9 1548  798 (S20W32 L=105)  1242    429        2 (1498/ 616-1628/ 963)  71.6 W27-W38 S31-S10
     10  764  824 (S10E18 L= 55)   236     49        0 ( 737/ 797- 792/ 843)   0.0 E19-E16 S11-S08
     11  233  716 (S10E62 L= 11)   808     79        1 ( 214/ 670- 259/ 755)  51.3 E65-E58 S12-S07
     12  938  698 (S19E09 L= 64)   199     89        0 ( 904/ 665- 979/ 723)  21.8 E11-E06 S20-S17
     13 1638  653 (S30W40 L=113)   681    143        1 (1574/ 618-1687/ 711)  26.6 W34-W44 S32-S26
     14  967  569 (S28E09 L= 63)   267     72        2 ( 946/ 528- 988/ 608) -90.0 E10-E08 S30-S25
     15  288  508 (S23E68 L=  5)   624     74        0 ( 270/ 477- 304/ 553)  26.6 E73-E63 S25-S20
     16  506  356 (S36E55 L= 18)   639    117        2 ( 462/ 314- 549/ 393)  13.2 E64-E47 S39-S34
  

    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