Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  743 1675 (N28E39 L=278)  1358    298        3 ( 662/1562- 830/1793) -24.0 E54-E27 N21-N34
      2 1076 1678 (N37E17 L=300)  1230    213        1 (1015/1583-1129/1761) -35.5 E25-E09 N31-N43
      3 1567 1670 (N56W26 L=343)  1269    147        1 (1510/1632-1634/1695)  28.4 W18-W42 N50-N63
      4  265  954 (S26E52 L=265)  1499    289        2 ( 193/ 832- 368/1028)  36.4 E62-E41 S35-S20
      5  660  696 (S32E14 L=303)   746    218        2 ( 562/ 677- 774/ 716)  21.8 E21-E06 S35-S28
      6 1500  578 (S16W43 L=  0)   346    103        0 (1448/ 556-1544/ 602)  45.0 W40-W45 S18-S13
      7 1178  528 (S28W25 L=342)   867    137        1 (1136/ 458-1212/ 599) -68.2 W22-W26 S33-S23
      8 1502  445 (S23W51 L=  9)   338    113        0 (1480/ 410-1527/ 493) -71.6 W49-W53 S25-S19
      9  416  447 (S59E41 L=276)   519     59        2 ( 387/ 429- 444/ 474)   0.0 E46-E36 S60-S56
     10 1320  415 (S31W40 L=357)   864    292        0 (1197/ 375-1427/ 452)  14.7 W31-W50 S33-S27
     11  535  373 (S62E24 L=294)  1940    229        2 ( 454/ 309- 627/ 449) -18.4 E37-E09 S65-S56
     12 1263  181 (S48W63 L= 20)   531     43        1 (1246/ 154-1277/ 209)  -7.6 W56-W71 S49-S45
  

    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