Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1685 (N58E18 L=  2)  1978    229        3 ( 593/1615- 786/1725)   0.0 E29-E04 N52-N61
      2 1233 1643 (N41W31 L= 51)  2066    464        3 (1026/1613-1439/1668) -22.1 W15-W47 N34-N47
      3 1619 1520 (N24W56 L= 76)   377     57        0 (1599/1490-1643/1550)  76.0 W54-W57 N21-N25
      4  250 1330 (N36E62 L=318)   768    114        0 ( 217/1276- 293/1366)  40.6 E65-E55 N32-N38
      5  471 1311 (N33E35 L=345)   319     63        0 ( 444/1290- 499/1333)  11.3 E37-E32 N31-N34
      6 1162 1222 (N16W12 L= 32)   875    243        1 (1090/1139-1210/1303) -68.2 W09-W13 N10-N22
      7 1800 1033 (S09W54 L= 74)   351     70        1 (1780/1000-1814/1074)  53.1 W52-W56 S10-S06
      8 1830  977 (S13W57 L= 77)   415    101        0 (1798/ 943-1853/1022)  26.6 W52-W60 S14-S10
      9 1587  685 (S25W32 L= 52)  1263    198        2 (1558/ 612-1637/ 773) -72.9 W30-W36 S31-S18
     10 1296  686 (S19W11 L= 31)   644    190        1 (1244/ 611-1369/ 740)  26.6 W06-W16 S22-S15
  

    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