Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  750 1598 (N19E34 L=134)  1144    167        4 ( 703/1543- 816/1662) -14.0 E40-E28 N15-N22
      2 1550 1560 (N43W19 L=188)   317     48        0 (1523/1539-1583/1582)   0.0 W16-W22 N42-N44
      3 1835 1154 (N27W50 L=218)   943    193        0 (1787/1084-1878/1223) -24.4 W44-W56 N24-N30
      4 1746  850 (N06W50 L=218)   333     86        0 (1710/ 833-1776/ 862)  15.9 W46-W53 N05-N07
      5 1073  800 (S17W09 L=177)   732    173        1 ( 983/ 754-1171/ 857)  55.0 W04-W13 S21-S10
      6 1762  778 (N03W55 L=223)   452     79        1 (1737/ 745-1778/ 812) -36.9 W52-W57 N01-N05
      7 1740  669 (S03W58 L=226)   342     45        1 (1720/ 646-1758/ 692)  90.0 W55-W59 S04-S01
      8 1264  606 (S22W27 L=195)   723    251        1 (1172/ 556-1367/ 677)   3.8 W19-W34 S22-S19
      9 1447  632 (S15W37 L=205)   232     58        1 (1436/ 602-1464/ 664) -45.0 W35-W39 S15-S12
     10 1714  580 (S08W60 L=228)   621     62        0 (1682/ 557-1738/ 605)  36.9 W57-W63 S09-S06
     11  788  499 (S43W02 L=170)  1135    307        1 ( 647/ 470- 946/ 537)  22.2 E09-W13 S47-S38
     12 1473  319 (S30W57 L=225)  1413    349        5 (1394/ 289-1575/ 361)  17.4 W48-W65 S32-S25
  

    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