Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1132 1728 (N46W34 L=158)   323     46        1 (1116/1692-1148/1758)  33.7 W31-W37 N43-N47
      2 1362 1563 (N28W39 L=162)   936    201        1 (1283/1504-1443/1602)   0.0 W30-W46 N26-N29
      3 1479 1414 (N17W39 L=162)   572    176        1 (1401/1377-1549/1458)  -4.4 W32-W45 N15-N17
      4 1060 1348 (N25W11 L=134)   758    162        2 (1003/1302-1136/1376) -12.5 W06-W15 N22-N26
      5 1772 1029 (S13W48 L=171)  5904    760       17 (1579/ 813-1887/1242)  18.4 W28-W67 S21-S03
      6 1029 1164 (N15W04 L=127)   226     85        0 ( 998/1123-1052/1200) -78.7 W02-W04 N12-N17
      7 1157 1099 (N08W09 L=132)   203    115        0 (1126/1056-1197/1146)  26.6 W06-W12 N06-N09
      8  841 1029 (N12E10 L=113)   189     55        0 ( 808/1015- 863/1048)   0.0 E12-E08 N10-N12
      9  543  969 (N15E30 L= 93)   472     85        1 ( 501/ 946- 583/1001) -33.7 E33-E27 N12-N16
     10  759  897 (N06E18 L=105)   302    101        1 ( 732/ 847- 775/ 935)  78.7 E19-E16 N03-N08
     11 1800  740 (S32W49 L=172)   785    163        0 (1739/ 703-1844/ 804)   8.7 W41-W54 S33-S29
     12  823  751 (S04E18 L=105)   512    162        0 ( 755/ 697- 897/ 793) -58.0 E21-E14 S08-N00
     13 1497  686 (S25W20 L=144)   681    126        3 (1463/ 642-1539/ 753) -83.7 W18-W23 S28-S19
     14 1506  553 (S34W19 L=142)  1752    392        1 (1385/ 429-1597/ 714)  31.8 W07-W28 S39-S25
     15 1153  554 (S23E04 L=119)   423    165        0 (1092/ 509-1207/ 615)  24.0 E09-W00 S24-S20
     16 1656  480 (S46W34 L=157)   982    116        1 (1610/ 439-1700/ 525)   5.2 W28-W39 S50-S42
     17  973  419 (S28E20 L=104)   417    123        0 ( 938/ 392-1025/ 471) -74.1 E21-E16 S30-S23
     18  422  398 (S18E62 L= 61)   569     50        1 ( 404/ 367- 437/ 430)  39.8 E65-E59 S20-S15
  

    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