Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  435 1501 (N22E50 L=181)  1001    299        2 ( 391/1431- 481/1587) -31.0 E58-E43 N17-N26
      2 1261 1453 (N31W10 L=241)   233     69        0 (1228/1417-1292/1475) -31.0 W07-W12 N28-N31
      3  968 1383 (N22E09 L=223)   202     65        0 ( 950/1352- 985/1412) -63.4 E09-E07 N19-N23
      4  591 1339 (N14E32 L=199)   341    160        0 ( 547/1293- 635/1383) -26.6 E36-E28 N11-N15
      5  961 1216 (N11E06 L=225)   288     61        1 ( 935/1187- 983/1255) -53.1 E08-E04 N09-N13
      6  323 1202 (N02E51 L=180)   296     61        0 ( 303/1172- 339/1224) -26.6 E53-E48 N00-N03
      7 1612  987 (N06W39 L=270)  2129    502        7 (1528/ 827-1720/1180)  67.8 W33-W46 S04-N18
      8 1501 1132 (N13W29 L=260)   499    197        0 (1397/1103-1568/1156)   0.0 W21-W34 N11-N14
      9  138  937 (S17E68 L=163)  2549    381        8 ( 104/ 831- 178/1022)  85.2 E74-E63 S23-S11
     10 1297  746 (S13W21 L=252)   973    309        2 (1182/ 670-1409/ 831)  42.5 W14-W27 S18-S06
     11 1654  674 (S13W49 L=280)  1235    155        2 (1618/ 602-1695/ 747)  84.3 W45-W51 S17-S07
     12 1362  686 (S16W26 L=258)   417    119        0 (1313/ 663-1431/ 703)  23.2 W23-W30 S17-S14
     13 1363  439 (S32W34 L=266)   756    142        2 (1288/ 422-1430/ 465)  21.8 W29-W39 S34-S29
  

    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