Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1691 (N38E24 L= 38)   428     76        1 ( 725/1667- 807/1720) -20.6 E28-E20 N36-N40
      2  374 1628 (N33E60 L=  3)  1506    165        1 ( 345/1548- 393/1687) -28.6 E76-E53 N26-N38
      3  466 1380 (N16E39 L= 23)  1372    402        2 ( 411/1221- 537/1556)  84.8 E42-E36 N05-N27
      4 1503 1449 (N25W31 L= 93)   333    117        1 (1449/1419-1554/1480) -14.0 W26-W34 N23-N26
      5  647 1327 (N13E25 L= 37)   746    254        0 ( 561/1257- 749/1387) -28.3 E31-E18 N08-N16
      6  566 1259 (N08E30 L= 33)   260    101        0 ( 532/1214- 580/1289) -59.0 E32-E28 N05-N10
      7 1938 1094 (N07W72 L=135)   956     67        4 (1916/1064-1964/1118) -26.6 W68-W79 N05-N08
      8  670  923 (S12E22 L= 41)   204     56        0 ( 642/ 908- 700/ 940)  14.0 E23-E19 S12-S10
      9 1245  881 (S12W14 L= 77)   441    205        1 (1182/ 812-1287/ 937)  58.0 W10-W16 S16-S08
     10  321  849 (S17E48 L= 14)  1389    304        4 ( 214/ 799- 476/ 882)   2.4 E59-E35 S19-S14
     11  562  779 (S21E30 L= 33)   223     45        1 ( 542/ 756- 581/ 812)  63.4 E31-E28 S22-S18
     12 1338  727 (S21W22 L= 84)   610    134        2 (1294/ 674-1381/ 779) -50.2 W18-W25 S23-S17
     13 1128  716 (S22W08 L= 71)   432    173        1 (1053/ 680-1201/ 750) -21.8 W03-W13 S24-S20
     14 1813  702 (S18W63 L=125)  1105    125        0 (1753/ 669-1870/ 736)  22.6 W56-W69 S20-S15
     15 1474  499 (S34W38 L=101)   645     75        2 (1434/ 471-1509/ 526) -20.6 W33-W41 S35-S32
     16  946  412 (S44E03 L= 60)   300     86        2 ( 911/ 389- 982/ 439) -21.8 E05-W00 S45-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