Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  430 1633 (N35E51 L=289)   424     56        1 ( 400/1609- 457/1661) -14.0 E54-E46 N33-N37
      2  621 1641 (N35E32 L=307)   334     41        0 ( 590/1623- 646/1656) -11.3 E34-E29 N33-N36
      3  949 1379 (N17E05 L=334)   258     82        0 ( 912/1352- 989/1401) -31.0 E07-E02 N15-N18
      4 1389 1242 (N10W22 L=  1)   571    186        1 (1338/1183-1461/1323) -53.1 W18-W26 N06-N14
      5  918 1153 (N04E07 L=333)   179     42        0 ( 902/1121- 939/1177)  76.0 E07-E05 N01-N05
      6 1603 1060 (N00W36 L= 16)   868    258        0 (1493/ 998-1704/1098)  14.0 W28-W44 S04-N02
      7 1814  876 (S09W56 L= 35)  1827    385        1 (1724/ 711-1860/1036)  78.1 W50-W60 S19-N00
      8  514  881 (S13E32 L=307)   223     62        1 ( 506/ 843- 531/ 920)  76.0 E33-E30 S14-S10
      9 1574  861 (S12W36 L= 15)   607    188        2 (1515/ 801-1607/ 900)  80.5 W31-W38 S15-S09
     10 1872  826 (S12W64 L= 43)  1783    166        4 (1826/ 764-1933/ 886)  28.6 W59-W72 S15-S08
     11  495  620 (S29E37 L=302)  2057    510        2 ( 382/ 528- 663/ 781) -21.8 E46-E24 S34-S18
     12 1557  753 (S18W36 L= 15)   275     63        2 (1523/ 734-1587/ 774)   0.0 W33-W38 S19-S16
     13  308  646 (S27E54 L=285)   576    155        0 ( 287/ 565- 327/ 717)  77.5 E56-E52 S31-S22
     14 1119  654 (S26W07 L=346)   854    204        1 (1016/ 626-1215/ 701)  13.0 W00-W13 S27-S22
     15 1614  550 (S31W46 L= 25)  1605    298        2 (1505/ 491-1723/ 655)  33.7 W37-W53 S34-S23
     16 1097  424 (S42W07 L=346)  2181    563        3 ( 799/ 396-1326/ 469)   5.7 E16-W25 S44-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