Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1525 1334 (N20W36 L=349)  2093    527        9 (1337/1226-1659/1407) -30.6 W23-W45 N12-N26
      2 1835 1286 (N14W65 L= 18)   954    125        3 (1801/1238-1864/1331) -78.7 W60-W67 N11-N16
      3  855 1225 (N18E10 L=303)   219    126        0 ( 813/1167- 888/1280) -54.5 E12-E07 N14-N21
      4  268 1084 (N10E54 L=259)   637    182        3 ( 223/1068- 326/1108)   5.2 E59-E48 N09-N11
      5  549 1053 (N08E31 L=282)   255     45        0 ( 522/1031- 569/1076)  18.4 E32-E29 N06-N09
      6 1873  940 (S07W65 L= 18)  1129    230        3 (1819/ 903-1915/ 968) -19.7 W57-W71 S10-S05
      7  646  909 (S01E24 L=289)   245     53        2 ( 616/ 885- 670/ 930) -33.7 E26-E22 S02-N00
      8 1530  750 (S16W32 L=345)  1198    330        1 (1414/ 634-1632/ 897)  45.0 W24-W39 S22-S07
      9  165  826 (S07E70 L=243)   815    114        2 ( 149/ 794- 194/ 872) -68.2 E72-E65 S08-S03
     10 1791  800 (S15W55 L=  8)   424    114        0 (1758/ 762-1809/ 823) -76.0 W51-W57 S17-S13
     11 1623  556 (S30W43 L=356)  1512    228        2 (1521/ 513-1698/ 605)  15.9 W34-W49 S32-S27
     12  448  533 (S25E46 L=268)   761    149        2 ( 403/ 504- 516/ 573) -26.6 E49-E40 S27-S22
  

    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