Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  750 1711 (N30E36 L= 66)  1194    233        1 ( 654/1643- 832/1784) -15.9 E48-E26 N26-N34
      2  426 1670 (N23E61 L= 40)  2486    224       11 ( 364/1565- 489/1722) -74.7 E72-E55 N15-N27
      3 1134 1640 (N33E07 L= 94)   774    166        2 (1089/1564-1183/1686) -33.7 E10-E01 N28-N36
      4 1332 1384 (N21W11 L=112)   851    318        4 (1221/1295-1442/1490) -22.4 W02-W19 N17-N25
      5 1102 1289 (N10E01 L=100)   942    462        2 (1064/1103-1136/1488) -64.5 E06-W04 N00-N21
      6 1776 1290 (N27W47 L=148)   832    107        1 (1751/1236-1802/1346) -63.4 W44-W49 N23-N29
      7 1331 1290 (N15W13 L=114)   200     75        0 (1312/1255-1344/1317) -63.4 W11-W13 N12-N16
      8  154  829 (S32E70 L= 31)   538     78        2 ( 119/ 811- 176/ 865)  12.1 E78-E64 S32-S29
      9 1326  814 (S13W22 L=123)   413     99        3 (1284/ 769-1359/ 857) -18.4 W18-W24 S14-S10
     10 1460  788 (S11W31 L=132)   515    189        2 (1388/ 753-1548/ 819)  31.0 W27-W37 S14-S07
     11 1221  765 (S18W17 L=118)   276    145        2 (1197/ 724-1266/ 818) -26.6 W14-W21 S19-S14
     12  966  723 (S25W03 L=104)   415    131        1 ( 916/ 674-1017/ 759) -14.0 E01-W07 S27-S23
     13 1715  648 (S12W56 L=157)   411     53        3 (1694/ 622-1741/ 671)   0.0 W53-W58 S13-S10
     14  871  611 (S34E01 L=100)   319    112        1 ( 826/ 587- 921/ 652) -15.9 E04-W03 S35-S32
     15 1487  418 (S31W51 L=152)   458    122        1 (1430/ 399-1541/ 446)  10.3 W44-W55 S31-S29
     16  837  200 (S68W15 L=116)   831    129        2 ( 772/ 167- 893/ 233)  -2.5 W03-W26 S70-S65
  

    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