Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1383 1552 (N36W07 L= 82)  2443    530        1 (1299/1323-1500/1765) -49.6 E05-W20 N24-N52
      2  373 1634 (N14E68 L=  7)  3239    157        7 ( 316/1574- 444/1696)  31.0 E80-E58 N09-N18
      3  432 1574 (N12E56 L= 19)   379     51        2 ( 408/1552- 454/1596)  90.0 E58-E53 N09-N13
      4  737 1512 (N15E31 L= 44)   486    166        1 ( 669/1471- 774/1541)  18.4 E33-E26 N11-N16
      5 1191 1477 (N25E03 L= 72)   391    147        0 (1136/1426-1243/1517) -14.0 E07-W01 N23-N26
      6  477 1435 (N04E45 L= 30)   385     72        0 ( 440/1419- 513/1455)  45.0 E47-E42 N02-N06
      7 1310 1320 (N19W08 L= 83)   239     91        0 (1279/1290-1345/1358) -21.8 W05-W11 N18-N20
      8 1627 1289 (N27W30 L=105)   283     59        0 (1605/1263-1651/1318) -26.6 W28-W32 N26-N28
      9  411 1264 (S07E44 L= 31)   567    131        2 ( 382/1192- 430/1315) -76.0 E46-E41 S11-S03
     10 1744 1271 (N30W41 L=116)   341     43        0 (1732/1241-1755/1303)  76.0 W39-W42 N28-N32
     11 1708 1194 (N24W39 L=113)   621    142        0 (1656/1164-1749/1227)  12.5 W33-W42 N21-N27
     12  973 1087 (S03E04 L= 70)   232     27        1 ( 951/1063- 991/1113) -56.3 E05-E02 S04-S01
     13  117 1039 (S25E73 L=  1)   755     57        2 (  92/1019- 151/1056)   8.7 E80-E66 S26-S24
     14  526  972 (S21E29 L= 46)   239     89        0 ( 487/ 950- 555/ 997) -11.3 E31-E26 S21-S19
     15  955  947 (S11E02 L= 73)   515    137        2 ( 894/ 904- 999/ 990)  60.3 E04-W00 S14-S07
     16 1603  774 (S02W41 L=116)  1370    265        2 (1573/ 645-1634/ 893) -47.7 W36-W47 S08-N04
     17 1780  820 (N06W54 L=129)   701    105        0 (1754/ 771-1807/ 875) -63.4 W52-W56 N02-N09
     18 1821  787 (N05W60 L=135)   454     56        1 (1807/ 757-1838/ 819) -76.0 W58-W62 N03-N07
  

    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