Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  912 1723 (N32E32 L=301)   799    161        2 ( 837/1686- 974/1751)  45.0 E34-E27 N27-N35
      2  738 1604 (N20E36 L=297)   379    103        2 ( 699/1556- 773/1648) -82.9 E37-E34 N15-N23
      3 1488 1567 (N42W14 L=347)   359     60        1 (1456/1548-1520/1593)  38.7 W11-W16 N39-N44
      4 1283 1099 (N05W12 L=345)   772    285        4 (1212/1029-1419/1191)  63.4 W08-W17 N00-N14
      5  596 1093 (S13E27 L=306)   657    248        0 ( 517/1044- 699/1132)  14.0 E31-E19 S16-S11
      6 1726 1008 (N15W44 L= 17)   722    189        0 (1668/ 953-1793/1086) -12.1 W37-W51 N13-N17
      7 1504 1014 (N08W28 L=  1)   241     46        0 (1490/ 988-1519/1049) -56.3 W26-W28 N06-N09
      8  248  863 (S35E54 L=279)   446     79        0 ( 215/ 846- 294/ 878)   5.2 E58-E47 S36-S34
      9  475  704 (S40E28 L=305)  1061    301        1 ( 358/ 666- 616/ 754)  25.6 E40-E17 S45-S33
     10 1524  499 (S19W50 L= 23)   646    142        1 (1463/ 476-1589/ 522)  40.6 W46-W54 S21-S15
     11 1293  437 (S30W38 L= 11)  1062    194        3 (1205/ 409-1389/ 474)  45.0 W32-W44 S34-S25
  

    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