Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1395 1755 (N58W27 L= 13)   489     74        1 (1366/1730-1419/1786)  30.3 W21-W33 N55-N62
      2  644 1749 (N43E50 L=296)   694    101        0 ( 593/1720- 682/1778)  51.3 E52-E46 N40-N46
      3  433 1631 (N31E61 L=285)  1406    130        4 ( 393/1573- 463/1697) -60.9 E68-E55 N27-N36
      4 1182 1599 (N40W03 L=349)   421    147        0 (1150/1550-1227/1658) -45.0 W00-W07 N37-N44
      5 1381 1433 (N31W20 L=  6)   477    156        0 (1345/1355-1421/1480) -49.4 W16-W22 N26-N33
      6  480 1302 (N11E40 L=306)   751    156        3 ( 461/1227- 508/1376) -83.7 E42-E36 N06-N15
      7  852 1274 (N14E14 L=332)   658    219        3 ( 809/1170- 883/1371) -63.4 E17-E10 N07-N19
      8 1841 1332 (N30W67 L= 53)   640     88        0 (1805/1299-1873/1362) -16.7 W61-W71 N27-N31
      9 1149 1145 (N09W06 L=352)   395    122        1 (1123/1084-1182/1197)  74.1 W05-W07 N05-N12
     10 1790 1068 (N12W55 L= 41)   656     87        3 (1746/1032-1818/1094) -20.6 W49-W57 N10-N13
     11  812 1025 (S02E13 L=333)   195     47        0 ( 787/1011- 840/1039)   0.0 E14-E11 S02-S01
     12 1067  905 (S06W04 L=350)   314     87        1 (1052/ 856-1084/ 945) -80.5 W02-W05 S09-S03
     13  667  780 (S19E20 L=326)   253     85        2 ( 645/ 744- 690/ 814)  90.0 E21-E18 S20-S16
     14  744  705 (S22E14 L=332)   271    106        2 ( 685/ 684- 788/ 722)   8.1 E18-E11 S23-S21
     15 1335  588 (S22W27 L= 13)  2552    582        6 (1120/ 488-1597/ 683)  31.8 W13-W43 S31-S13
     16  556  572 (S34E29 L=317)   339    105        0 ( 539/ 524- 577/ 618) -80.5 E30-E27 S37-S31
     17  787  555 (S32E10 L=336)   525    136        1 ( 726/ 536- 857/ 580)   0.0 E14-E05 S33-S30
  

    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