Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1531 1325 (N20W37 L=350)  1862    390        4 (1359/1224-1653/1414) -34.4 W25-W44 N12-N26
      2 1833 1284 (N14W64 L= 18)   921    118        3 (1799/1236-1864/1333) -90.0 W60-W68 N10-N17
      3  842 1226 (N18E11 L=303)   193    117        0 ( 797/1170- 875/1278) -54.5 E13-E08 N14-N21
      4  885 1123 (N11E08 L=305)   240     71        1 ( 857/1101- 907/1157)  33.7 E10-E06 N10-N13
      5  264 1085 (N10E55 L=259)   666    177        2 ( 215/1069- 321/1109)   4.8 E60-E48 N09-N11
      6  537 1054 (N08E31 L=282)   314     44        1 ( 506/1033- 561/1075)  26.6 E33-E29 N06-N09
      7 1867  940 (S07W64 L= 17)  1114    231        2 (1812/ 902-1911/ 969) -24.8 W57-W71 S10-S04
      8  641  910 (S01E25 L=289)   248     47        1 ( 618/ 885- 666/ 932) -33.7 E26-E22 S02-N00
      9 1522  750 (S16W32 L=345)  1195    317        2 (1395/ 630-1623/ 893)  43.2 W23-W39 S22-S07
     10  162  825 (S07E70 L=243)   684    117        1 ( 145/ 802- 190/ 860) -31.0 E73-E65 S08-S04
     11 1616  556 (S30W42 L=356)  1424    264        3 (1514/ 518-1694/ 600)  15.9 W33-W49 S32-S27
     12  427  542 (S25E47 L=267)   788    140        2 ( 372/ 504- 481/ 581) -26.6 E51-E43 S27-S21
  

    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