Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1438 1539 (N29W42 L=229)   693    267        0 (1330/1498-1561/1583) -35.5 W35-W49 N23-N33
      2 1639 1506 (N21W57 L=244)   586     75        2 (1611/1478-1665/1537)  15.9 W53-W60 N20-N23
      3 1068 1356 (N26W09 L=197)  1190    249        3 ( 981/1260-1148/1448) -71.6 W05-W12 N18-N33
      4  923 1412 (N32W01 L=189)   234    120        0 ( 880/1379- 970/1437)  18.4 E02-W04 N30-N33
      5 1281 1380 (N23W24 L=211)   397    173        0 (1201/1363-1356/1400) -10.3 W18-W29 N21-N24
      6  649 1105 (N18E22 L=165)   281     80        0 ( 623/1066- 671/1145) -68.2 E23-E20 N15-N20
      7  876 1060 (N11E08 L=179)   278    179        2 ( 836/1012- 962/1088)  12.5 E11-E02 N09-N13
      8  993  824 (S04E05 L=182)   488    103        1 ( 947/ 791-1045/ 857) -50.2 E07-E02 S07-S01
      9 1123  678 (S16E00 L=187)   449    155        0 (1082/ 617-1145/ 750)  60.3 E04-W02 S18-S11
     10 1380  595 (S27W15 L=202)  2331    529        4 (1179/ 476-1566/ 749)   7.9 W00-W29 S35-S18
     11  355  523 (S14E60 L=127)   778    107        0 ( 317/ 475- 388/ 553) -59.0 E63-E57 S17-S12
     12  858  487 (S24E22 L=166)   287     58        0 ( 827/ 458- 883/ 511)  21.8 E24-E19 S25-S22
  

    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