Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1237 1811 (N54E15 L=244)   525     76        1 (1196/1782-1270/1831)  26.6 E18-E11 N51-N56
      2  855 1795 (N38E41 L=219)  1188    173        3 ( 777/1772- 938/1822)  28.6 E46-E33 N33-N40
      3  501 1704 (N22E61 L=198)   675    120        1 ( 476/1669- 538/1753) -66.8 E64-E57 N19-N26
      4  311 1549 (N09E69 L=190)  1020    112        4 ( 283/1499- 348/1579)  71.6 E73-E65 N05-N11
      5 1311 1356 (N22W08 L=267)   539    191        0 (1265/1320-1385/1409)  12.5 W03-W12 N19-N25
      6  219 1290 (S08E65 L=195)   676    115        2 ( 173/1273- 274/1312)   8.1 E71-E57 S10-S06
      7  272  957 (S27E51 L=209)   311     42        2 ( 252/ 935- 295/ 980) -14.0 E53-E48 S28-S25
      8 1843  875 (N11W60 L=320)  1324    164        4 (1792/ 811-1868/ 959)  68.2 W56-W65 N06-N16
      9 1501  597 (S15W42 L=302)   296     73        0 (1474/ 567-1526/ 634) -76.0 W41-W43 S16-S12
     10 1363  538 (S22W36 L=295)   464    242        0 (1314/ 504-1442/ 627)  -4.8 W29-W41 S23-S18
     11  732  562 (S40E05 L=254)   330    112        0 ( 710/ 516- 771/ 608) -29.7 E07-W00 S41-S37
  

    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