Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1860 (N52E22 L= 95)   612     36        2 (1093/1830-1138/1884) -39.8 E26-E19 N48-N54
      2  509 1684 (N22E55 L= 63)   507     52        0 ( 481/1668- 542/1703)  26.6 E57-E52 N20-N23
      3  382 1525 (N10E56 L= 61)   493     73        1 ( 345/1492- 415/1553)  59.0 E58-E53 N07-N12
      4 1732 1260 (N25W41 L=159)   618    129        0 (1704/1194-1747/1316) -69.4 W38-W42 N20-N28
      5  448  672 (S41E34 L= 84)   313     44        1 ( 423/ 653- 470/ 696)  45.0 E36-E31 S42-S38
      6 1713  648 (S10W56 L=173)   356     37        2 (1690/ 627-1735/ 671)  45.0 W53-W58 S11-S08
      7 1358  503 (S28W37 L=154)   547    127        0 (1314/ 459-1416/ 557)  76.0 W34-W38 S31-S23
      8  901  249 (S60W21 L=139)   448     58        1 ( 871/ 218- 930/ 283) -15.3 W15-W27 S62-S57
      9 1075  188 (S58W45 L=162)   732     63        2 (1049/ 159-1099/ 224) -15.9 W38-W52 S60-S55
  

    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