Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1500 1694 (N47W31 L=259)  1406    146        3 (1431/1652-1562/1733)   0.0 W23-W37 N43-N50
      2  652 1587 (N25E33 L=194)   630    142        0 ( 631/1522- 668/1649) -56.3 E37-E30 N20-N29
      3  391 1537 (N20E54 L=174)   932    161        2 ( 365/1473- 417/1614) -51.3 E58-E50 N15-N25
      4 1367 1393 (N21W17 L=245)  1471    397        1 (1193/1274-1490/1485) -20.9 W05-W26 N15-N25
      5  703 1427 (N15E25 L=202)   342     55        1 ( 674/1400- 733/1445) -14.0 E27-E23 N13-N16
      6 1823 1306 (N24W57 L=284)   375     50        2 (1796/1281-1840/1330)  56.3 W53-W59 N22-N25
      7 1762 1273 (N20W49 L=277)   510     59        1 (1740/1242-1783/1306) -76.0 W46-W51 N18-N22
      8 1858 1266 (N22W61 L=288)   551     57        3 (1835/1244-1883/1292) -45.0 W57-W64 N20-N23
      9  175 1259 (N02E67 L=160)   613     62        1 ( 157/1236- 200/1291) -20.6 E71-E63 N00-N03
     10  490 1038 (S10E34 L=194)   213     60        0 ( 465/1005- 513/1061)  53.1 E35-E32 S12-S08
     11  894  855 (S17E06 L=222)   376    187        0 ( 816/ 832- 987/ 881)  21.8 E10-W00 S19-S15
     12  312  844 (S23E49 L=178)   549    120        2 ( 263/ 825- 366/ 865)   0.0 E54-E43 S24-S21
     13  638  685 (S31E22 L=205)   463     91        1 ( 602/ 646- 674/ 725)  45.0 E25-E19 S33-S27
     14 1501  549 (S27W41 L=269)  2459    501        2 (1310/ 439-1692/ 690)  47.7 W30-W51 S37-S15
     15 1787  574 (S20W69 L=297)   565     67        1 (1761/ 550-1809/ 603)  45.0 W66-W72 S21-S18
     16 1425  357 (S41W47 L=274)   333     63        0 (1397/ 333-1446/ 386) -78.7 W45-W48 S43-S38
  

    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