Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1450 1764 (N46W54 L=127)  2064    268        1 (1402/1677-1510/1839)  61.7 W50-W60 N38-N51
      2 1795 1359 (N13W63 L=136)  1767    316        1 (1716/1228-1877/1469) -62.1 W58-W69 N03-N20
      3 1734 1353 (N14W56 L=129)   378     65        2 (1709/1327-1754/1378)  45.0 W53-W58 N12-N15
      4  435 1337 (N32E39 L= 34)   358    116        1 ( 391/1321- 480/1359) -18.4 E43-E34 N30-N33
      5  246 1284 (N30E57 L= 17)  2317    477        4 ( 136/1255- 344/1324)   5.0 E69-E46 N27-N32
      6 1472 1259 (N13W32 L=105)   284    116        0 (1423/1222-1509/1298)  23.2 W28-W35 N11-N15
      7  542 1069 (N14E31 L= 43)   338     93        0 ( 496/1038- 581/1102) -33.7 E33-E27 N11-N16
      8 1718  911 (S12W45 L=119)   366     55        1 (1697/ 881-1741/ 932) -76.0 W43-W47 S13-S09
      9 1786  675 (S29W53 L=126)  1444    218        3 (1712/ 610-1840/ 748)  28.6 W45-W58 S32-S24
     10 1666  581 (S33W42 L=115)   397     74        0 (1639/ 563-1693/ 602) -38.7 W38-W44 S34-S30
     11 1601  557 (S33W36 L=109)   440     67        2 (1571/ 536-1628/ 581)  36.9 W33-W38 S34-S31
     12  937  505 (S26E13 L= 60)  1000    242        1 ( 838/ 471-1033/ 566)   0.0 E19-E05 S28-S22
  

    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