Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1611 (N44E06 L=277)   304     76        0 ( 914/1580- 948/1649) -78.7 E06-E03 N41-N46
      2 1114 1593 (N41W10 L=293)   270     84        0 (1088/1570-1150/1623) -45.0 W07-W12 N39-N43
      3  976 1438 (N31E02 L=281)   248     43        2 ( 958/1413-1000/1464)  71.6 E03-W00 N29-N32
      4 1722 1276 (N16W52 L=335)   604    126        1 (1667/1235-1773/1307) -32.0 W47-W56 N13-N18
      5 1456 1196 (N13W29 L=312)   541    156        3 (1390/1148-1511/1250) -45.0 W24-W32 N09-N16
      6  231  848 (S06E60 L=223)   700     96        4 ( 195/ 823- 283/ 882) -18.4 E64-E54 S07-S04
      7  615  853 (S05E27 L=256)   436    147        1 ( 540/ 814- 668/ 880) -20.6 E31-E23 S07-S03
      8  350  716 (S14E50 L=233)   325     78        1 ( 330/ 685- 375/ 750)  63.4 E52-E46 S16-S12
      9  665  603 (S22E26 L=257)   333    165        1 ( 624/ 551- 691/ 658) -82.9 E28-E24 S25-S17
     10  292  588 (S23E62 L=221)   675    157        3 ( 258/ 525- 320/ 644)  60.9 E67-E58 S28-S19
     11  859  477 (S31E14 L=269)   258     89        2 ( 828/ 446- 897/ 508)  45.0 E16-E10 S33-S28
     12 1267  412 (S38W17 L=300)   589    135        2 (1239/ 348-1301/ 474)  73.3 W15-W18 S43-S33
  

    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