Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1495 1715 (N30W58 L=127)  1566    242        2 (1417/1620-1555/1789)  61.7 W52-W63 N23-N36
      2  673 1754 (N62E03 L= 66)   547     75        0 ( 645/1712- 694/1781) -90.0 E06-W00 N57-N65
      3  379 1641 (N62E49 L= 20)   813     72        1 ( 339/1615- 418/1665)   3.6 E61-E41 N60-N64
      4 1586 1618 (N22W57 L=126)   375     36        1 (1560/1599-1610/1634) -11.3 W55-W60 N20-N22
      5 1303 1515 (N24W31 L=100)   321     79        0 (1259/1499-1345/1528) -21.8 W28-W33 N22-N25
      6 1098 1425 (N25W15 L= 84)  1127    227        3 ( 986/1383-1212/1447) -24.8 W09-W22 N21-N28
      7 1375 1193 (N05W24 L= 93)   269     89        0 (1334/1173-1420/1220) -45.0 W21-W26 N03-N07
      8 1631 1127 (S05W38 L=107)   339     78        1 (1612/1088-1650/1170)  45.0 W36-W40 S07-S03
      9 1866 1036 (S18W58 L=127)   712    114        2 (1844/ 986-1904/1099)  29.1 W54-W64 S20-S15
     10  686 1046 (N14E19 L= 50)   293     86        0 ( 640/1025- 724/1073)   0.0 E21-E15 N12-N14
     11  941  982 (N04E06 L= 63)   255     25        0 ( 917/ 961- 963/1005)  14.0 E07-E03 N03-N05
     12  722  735 (S04E24 L= 45)  1165    339        1 ( 574/ 647- 825/ 847)  14.0 E35-E15 S06-N00
     13  240  715 (N04E61 L=  8)   710    121        2 ( 194/ 681- 299/ 756) -33.7 E65-E56 N01-N07
     14  372  497 (S11E59 L= 10)   629    126        1 ( 351/ 442- 395/ 548)  36.9 E63-E54 S13-S07
     15 1509  481 (S42W18 L= 87)  1099    255        2 (1403/ 401-1605/ 548)  11.3 W07-W27 S44-S38
     16  456  452 (S15E54 L= 15)   321     58        0 ( 439/ 426- 478/ 487)  63.4 E55-E52 S16-S12
     17  422  400 (S18E61 L=  8)   555     46        1 ( 393/ 380- 448/ 425) -21.8 E64-E58 S18-S15
     18 1256  339 (S44E06 L= 63)   363     90        0 (1211/ 322-1308/ 354) -15.9 E08-E01 S45-S41
     19 1020  251 (S43E27 L= 42)   670    127        1 ( 956/ 234-1083/ 275) -11.3 E31-E21 S44-S40
  

    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