Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  557 1644 (N38E47 L=328)   374     53        1 ( 537/1613- 581/1666) -33.7 E51-E44 N35-N39
      2  706 1625 (N38E32 L=343)  1061    168        3 ( 629/1586- 785/1662)   0.0 E39-E24 N36-N40
      3 1278 1391 (N27W14 L= 30)  1314    282        2 (1223/1294-1331/1500) -81.9 W10-W18 N20-N34
      4  898 1470 (N29E13 L=  3)   249     48        0 ( 872/1453- 922/1494)   0.0 E14-E10 N27-N30
      5  451 1364 (N18E44 L=332)   535    168        0 ( 400/1299- 486/1411) -28.6 E49-E38 N14-N20
      6 1584 1283 (N22W37 L= 53)   483    133        3 (1548/1236-1631/1340)  36.9 W34-W42 N18-N25
      7  902 1119 (N07E08 L=  7)   760    258        4 ( 819/1049-1013/1184)  42.3 E12-E01 N01-N11
      8  679  863 (S11E21 L=355)   205     64        0 ( 661/ 846- 712/ 885) -33.7 E21-E18 S11-S09
      9 1270  565 (S25W21 L= 36)   410     84        3 (1226/ 539-1317/ 593)  -8.1 W17-W24 S26-S23
     10 1588  512 (S27W48 L= 64)   347     85        0 (1551/ 491-1631/ 541)  33.7 W45-W51 S28-S24
  

    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