Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1670 1546 (N41W51 L=127)   990    157        2 (1610/1483-1727/1593) -45.0 W45-W54 N36-N44
      2 1075 1492 (N32E01 L= 76)   341     51        1 (1049/1473-1105/1515) -33.7 E02-W01 N30-N33
      3 1211 1412 (N28W10 L= 86)   406     71        1 (1189/1379-1246/1448) -78.7 W08-W12 N25-N30
      4  283 1287 (N11E57 L= 20)  1961    316        2 ( 254/1120- 319/1401) -76.8 E61-E52 N01-N18
      5  152 1258 (N08E75 L=  2)  1063    123        2 ( 130/1212- 171/1307)  90.0 E78-E71 N05-N11
      6  974  840 (S10E02 L= 75)   553    180        1 ( 954/ 757-1006/ 931) -78.7 E03-W00 S14-S04
      7 1750  722 (S12W56 L=132)   598    123        2 (1720/ 685-1794/ 747)  38.7 W52-W59 S14-S10
      8 1325  652 (S19W23 L= 99)   326     94        0 (1306/ 613-1352/ 691) -68.2 W21-W25 S21-S16
      9 1456  568 (S24W34 L=111)   260     44        0 (1437/ 544-1476/ 594)  90.0 W32-W35 S25-S21
     10 1369  463 (S32W31 L=107)   423     96        1 (1324/ 448-1417/ 478)  23.2 W27-W34 S33-S30
  

    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