Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1192 1843 (N58E19 L=218)   981     66        3 (1164/1804-1217/1879) -63.4 E25-E15 N53-N63
      2 1392 1604 (N43W10 L=247)   404     90        1 (1365/1563-1428/1636) -23.2 W06-W13 N40-N45
      3 1595 1602 (N51W31 L=268)   367     35        1 (1576/1575-1615/1621) -56.3 W28-W33 N48-N53
      4 1763 1477 (N47W52 L=290)   670     68        1 (1738/1441-1797/1501)   0.0 W47-W58 N43-N48
      5 1411 1335 (N25W16 L=254)   284    122        0 (1380/1294-1448/1391) -33.7 W12-W19 N23-N27
      6 1764 1183 (N25W47 L=285)  2035    385        1 (1662/1031-1859/1316) -30.6 W36-W59 N18-N31
      7  265  952 (S25E52 L=186)   993    196        1 ( 222/ 869- 313/1028)  74.7 E57-E47 S30-S19
      8 1235  670 (S18W22 L=260)   236    100        0 (1187/ 656-1289/ 683)  26.6 W18-W24 S19-S16
      9 1181  515 (S29W25 L=262)  2037    469        4 ( 969/ 444-1383/ 581)  39.3 W13-W35 S38-S19
     10 1464  482 (S23W45 L=283)   376     63        0 (1438/ 458-1494/ 511)  76.0 W44-W46 S24-S20
     11 1495  306 (S32W64 L=302)  1356    135        1 (1430/ 288-1568/ 330)   8.1 W56-W77 S34-S30
     12 1034  151 (S61W50 L=288)   690     80        0 (1017/ 120-1056/ 183) -10.8 W42-W63 S63-S57
  

    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