Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1144 1612 (N32E00 L= 47)   250     70        0 (1107/1596-1184/1627) -14.0 E02-W02 N31-N33
      2  312 1453 (N13E56 L=351)  4251    307       13 ( 233/1338- 411/1550)  60.3 E64-E48 N05-N19
      3 1722 1501 (N35W47 L= 95)   551     82        0 (1684/1474-1751/1524)   7.1 W42-W50 N33-N36
      4  543 1471 (N16E37 L= 10)   339     78        2 ( 510/1441- 569/1505)  71.6 E39-E35 N13-N17
      5  909 1131 (S01E08 L= 39)   197     61        0 ( 890/1099- 925/1163) -63.4 E09-E06 S02-N00
      6  253 1039 (S12E53 L=354)  1480    254        3 ( 204/ 947- 306/1152)  67.4 E58-E48 S17-S05
      7 1859 1075 (N10W58 L=105)   943    124        2 (1838/1009-1878/1132)  69.4 W55-W60 N05-N13
      8  166 1080 (S10E63 L=344)  1094    117        2 ( 132/1031- 189/1125)  90.0 E68-E60 S12-S07
      9 1539  755 (S14W36 L= 83)  1139    452        1 (1463/ 610-1629/ 924)  90.0 W32-W41 S23-S03
     10 1004  837 (S17W01 L= 48)   559    174        0 ( 974/ 752-1029/ 917) -73.3 E01-W03 S22-S12
     11  278  846 (S24E51 L=356)  3163    538        8 (  99/ 771- 572/ 900)   7.1 E82-E26 S28-S18
     12 1236  848 (S13W15 L= 62)   289    127        0 (1178/ 831-1299/ 867)   0.0 W11-W19 S14-S12
     13  587  686 (S31E26 L= 21)  2592    522        4 ( 389/ 612- 756/ 859)  31.8 E43-E14 S36-S18
     14  952  426 (S44W04 L= 51)   331     32        0 ( 933/ 402- 971/ 454) -56.3 W02-W06 S46-S42
     15 1161  401 (S42W21 L= 69)  1981    351        4 ( 980/ 366-1328/ 437)   7.9 W06-W35 S44-S39
     16 1601  369 (S36W63 L=110)   728     78        2 (1565/ 350-1641/ 394)  12.5 W57-W67 S36-S33
  

    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