Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1754 (N58W02 L= 62)   370     33        1 ( 853/1731- 905/1771)   0.0 W00-W05 N56-N60
      2  681 1720 (N59E22 L= 38)   405     43        1 ( 649/1702- 704/1739)  15.9 E26-E18 N57-N60
      3 1057 1438 (N31W09 L= 68)   361    131        0 (1003/1389-1085/1463) -56.3 W05-W10 N27-N33
      4  850 1365 (N29E07 L= 53)   417    122        1 ( 801/1327- 907/1410) -39.8 E09-E02 N26-N32
      5 1552 1306 (N15W39 L= 99)   251     59        1 (1528/1277-1579/1325) -26.6 W37-W41 N12-N15
      6 1560 1240 (N11W38 L= 97)   339     59        0 (1525/1217-1590/1262) -26.6 W35-W39 N08-N12
      7  408 1165 (N22E41 L= 19)   634    189        1 ( 388/1087- 424/1254) -74.7 E42-E39 N16-N27
      8 1915  971 (S13W69 L=128)  2126    192        2 (1886/ 898-1947/1054) -78.7 W64-W77 S18-S08
      9 1556  959 (S06W32 L= 92)   272     74        0 (1544/ 919-1569/1002) -90.0 W31-W33 S08-S03
     10 1220  748 (S13W08 L= 68)   344    112        0 (1180/ 726-1272/ 784)  18.4 W05-W11 S14-S11
     11 1769  629 (S32W53 L=113)  1498    204        0 (1697/ 562-1858/ 686)   8.1 W45-W67 S35-S29
     12  389  584 (S14E52 L=  8)   799    191        1 ( 344/ 506- 414/ 665)  42.5 E60-E47 S19-S08
     13  245  602 (S12E68 L=351)   592     62        1 ( 221/ 574- 270/ 619) -31.0 E71-E63 S14-S11
     14 1278  364 (S42W08 L= 68)   511    103        0 (1238/ 341-1333/ 381)   0.0 W04-W13 S43-S39
     15 1091  269 (S47E09 L= 50)   996    122        1 (1039/ 232-1159/ 294)   9.5 E14-E02 S49-S43
     16  920  204 (S50E28 L= 32)  1303    133        1 ( 880/ 156- 960/ 257)  38.2 E36-E21 S55-S44
  

    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