Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  841 1515 (N29E20 L= 51)  1548    318        5 ( 728/1385- 956/1632) -30.6 E31-E09 N21-N35
      2 1172 1598 (N40W02 L= 73)   241     35        0 (1150/1581-1197/1614)   0.0 W00-W04 N38-N41
      3  521 1421 (N18E40 L= 31)  2173    451        6 ( 436/1278- 578/1563) -41.6 E51-E32 N10-N26
      4 1224 1474 (N32W08 L= 79)   292     99        0 (1209/1428-1238/1527) -74.1 W06-W08 N28-N35
      5 1026 1330 (N19E04 L= 67)   215    112        0 ( 970/1312-1084/1347)  15.9 E07-W00 N17-N21
      6 1284 1243 (N17W14 L= 85)   892    307        2 (1190/1182-1348/1318)  16.7 W07-W17 N13-N21
      7 1512 1220 (N19W30 L=101)   274     90        0 (1499/1167-1528/1266) -71.6 W28-W30 N15-N21
      8 1078 1186 (N11W01 L= 72)   274     66        0 (1044/1157-1112/1219) -26.6 E01-W03 N09-N12
      9 1763 1086 (N13W51 L=122)   702    131        1 (1741/1031-1792/1138)  66.8 W49-W53 N09-N16
     10  859 1126 (N05E11 L= 60)   265    116        0 ( 802/1112- 920/1136)  15.9 E14-E07 N03-N05
     11  206  982 (S12E60 L= 11)   862    178        3 ( 160/ 955- 257/1018) -15.3 E66-E53 S13-S10
     12 1324  401 (S36W33 L=104)  1757    329        6 (1194/ 373-1475/ 437)  12.8 W22-W45 S39-S31
  

    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