Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1646 (N38E48 L=328)   368     49        1 ( 530/1617- 572/1666) -23.2 E52-E45 N36-N39
      2  700 1627 (N38E33 L=343)  1127    185        3 ( 622/1589- 783/1665)  -3.6 E40-E24 N36-N41
      3  303 1520 (N27E69 L=307)   648     58        1 ( 273/1496- 337/1540)  -7.6 E77-E62 N25-N27
      4 1273 1393 (N27W14 L= 30)  1316    298        1 (1217/1293-1324/1497) -81.3 W10-W17 N20-N34
      5  893 1470 (N29E13 L=  3)   247     35        1 ( 869/1453- 915/1492)  18.4 E14-E11 N27-N30
      6  447 1364 (N18E44 L=332)   587    156        0 ( 393/1302- 478/1408) -31.0 E50-E39 N14-N20
      7 1581 1284 (N22W37 L= 53)   382    130        1 (1542/1245-1631/1339)  36.9 W33-W42 N19-N25
      8 1531 1200 (N16W33 L= 48)   243     41        0 (1509/1174-1551/1221) -56.3 W30-W34 N14-N17
      9  895 1123 (N07E09 L=  7)   901    258        4 ( 817/1053-1008/1188)  42.0 E12-E02 N02-N11
     10  695  855 (S11E20 L=356)   338    122        0 ( 656/ 823- 764/ 888) -20.6 E22-E14 S12-S09
     11 1265  564 (S25W21 L= 36)   389     74        1 (1222/ 543-1305/ 591) -18.4 W17-W23 S26-S23
     12 1124  387 (S40W14 L= 30)   290     57        1 (1091/ 373-1153/ 413)  21.8 W11-W16 S41-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