Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1655 (N47W14 L=161)   363    108        0 ( 904/1614- 977/1694)  18.4 W10-W19 N45-N49
      2 1035 1398 (N28W11 L=157)   450    123        1 ( 974/1359-1082/1427)   6.3 W05-W14 N26-N29
      3 1158 1248 (N16W13 L=160)   214    122        0 (1099/1232-1207/1271) -18.4 W10-W16 N15-N18
      4 1894 1172 (S09W67 L=213)  2800    150        4 (1843/1110-1938/1259)  14.0 W58-W75 S13-S04
      5  717 1203 (N24E14 L=132)   826    206        2 ( 630/1145- 803/1258)   0.0 E21-E08 N22-N26
      6  423 1106 (N25E37 L=109)   430    135        1 ( 369/1063- 477/1143)  11.3 E42-E32 N22-N26
      7 1513  846 (S15W25 L=171)  2382    619       10 (1318/ 717-1649/1082)  22.5 W09-W39 S19-S05
      8 1835  948 (S20W53 L=200)   653    152        0 (1801/ 890-1857/1027)  36.9 W49-W57 S22-S15
      9  275  896 (N15E54 L= 93)   300     53        0 ( 262/ 867- 288/ 930)  76.0 E55-E52 N12-N16
     10 1525  911 (S12W27 L=173)   310     51        0 (1497/ 886-1550/ 929) -33.7 W24-W28 S13-S10
     11  540  609 (S06E40 L=107)   854    245        0 ( 439/ 551- 635/ 688) -50.7 E45-E34 S11-N00
     12 1544  603 (S32W23 L=170)   517    121        0 (1484/ 571-1591/ 635)   0.0 W18-W27 S32-S30
     13 1195  540 (S25E02 L=145)   377     74        1 (1152/ 508-1221/ 566) -36.9 E04-W00 S27-S23
     14 1103  485 (S27E09 L=138)   374    148        1 (1073/ 431-1152/ 538)  32.0 E12-E04 S29-S23
     15  828  435 (S23E28 L=119)   457    133        0 ( 791/ 392- 863/ 499)  32.0 E31-E23 S25-S19
     16 1170  392 (S35E08 L=139)   918    195        3 (1095/ 323-1231/ 473)  31.6 E15-E02 S38-S30
     17  663  373 (S24E42 L=104)  1008    221        3 ( 614/ 303- 732/ 457) -90.0 E46-E38 S29-S17
     18 1426  368 (S46W10 L=157)   551     90        1 (1381/ 342-1461/ 399)   0.0 W06-W13 S47-S43
     19 1061  137 (S56E35 L=111)   519     52        0 (1042/ 114-1078/ 169)  32.5 E40-E29 S59-S52
  

    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