Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  668 1828 (N40E59 L= 60)  1355    105        0 ( 633/1773- 704/1881) -12.7 E79-E48 N35-N43
      2  886 1664 (N33E27 L= 92)   408     79        0 ( 863/1624- 914/1703) -74.1 E28-E25 N29-N36
      3  334 1548 (N13E65 L= 54)  3005    237        6 ( 261/1487- 386/1664) -81.3 E73-E56 N08-N21
      4  425 1601 (N18E59 L= 60)   757     59        1 ( 398/1570- 451/1628) -20.6 E62-E54 N16-N20
      5  247 1351 (N00E63 L= 56)   434     87        0 ( 237/1308- 259/1390) -38.7 E66-E60 S02-N02
      6  992 1095 (N01E03 L=116)  1360    343        1 ( 949/ 935-1044/1216) -52.4 E08-W03 S07-N07
      7 1794 1082 (N19W51 L=170)  1817    284        1 (1708/ 969-1864/1188) -21.8 W41-W61 N13-N23
      8  475 1111 (S09E36 L= 83)   492    140        0 ( 412/1094- 538/1142)  24.0 E40-E31 S11-S06
      9 1122  716 (S18W13 L=132)   223     49        0 (1107/ 686-1141/ 746) -45.0 W11-W14 S19-S16
     10 1573  606 (S14W46 L=165)  1763    349        0 (1414/ 546-1720/ 669)  37.9 W36-W56 S20-S06
     11 1580  411 (S24W59 L=178)  1996    183        3 (1488/ 385-1674/ 439)  14.7 W50-W70 S27-S20
     12  930  377 (S45W13 L=132)   593     78        1 ( 890/ 354- 975/ 409)   8.1 W09-W16 S46-S42
  

    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