Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1487 (N32E28 L=190)   689    158        1 ( 497/1460- 632/1521) -15.3 E33-E22 N29-N34
      2 1565 1478 (N22W41 L=260)   277     74        0 (1524/1465-1593/1494) -11.3 W38-W43 N21-N22
      3 1231 1372 (N18W16 L=234)   345     85        1 (1179/1348-1282/1401) -18.4 W12-W18 N15-N19
      4  435 1337 (N23E37 L=181)   576    153        0 ( 388/1279- 467/1379) -48.8 E41-E34 N18-N26
      5 1455 1159 (N03W28 L=246)   390     83        0 (1439/1108-1470/1207)  80.5 W26-W28 N00-N06
      6 1718 1108 (S02W47 L=265)   856    146        1 (1676/1074-1775/1162) -26.6 W43-W51 S04-N01
      7  495 1087 (N07E33 L=186)   410     47        1 ( 468/1056- 516/1117)  76.0 E34-E31 N04-N08
      8 1745  972 (S10W48 L=266)   396     46        0 (1724/ 939-1764/1005)  90.0 W46-W49 S11-S07
      9  333  962 (N01E46 L=172)   286     72        1 ( 300/ 943- 351/ 990)  21.8 E49-E44 N00-N02
     10 1548  780 (S20W32 L=251)  1975    592        4 (1418/ 547-1679/ 938)  49.4 W23-W42 S34-S10
     11  261  816 (S07E55 L=163)  2102    366        2 ( 196/ 697- 343/ 894) -63.4 E62-E49 S14-S01
     12  650  837 (S09E25 L=194)   265     47        1 ( 628/ 806- 665/ 865) -90.0 E26-E23 S11-S07
     13 1445  740 (S22W25 L=243)   378    176        0 (1372/ 705-1510/ 805)  24.0 W19-W29 S23-S18
     14 1329  765 (S19W17 L=235)   192     51        1 (1296/ 745-1359/ 786)  11.3 W14-W19 S20-S18
     15  383  622 (S20E49 L=170)   689    159        3 ( 360/ 559- 407/ 682)  45.0 E52-E45 S23-S16
     16  632  558 (S26E31 L=187)   492     55        2 ( 598/ 531- 664/ 581) -26.6 E33-E28 S28-S24
     17 1561  412 (S46W40 L=258)   881     92        3 (1532/ 373-1584/ 456) -54.5 W37-W43 S49-S42
  

    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