Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  541 1709 (N47E49 L=130)  3400    462        7 ( 425/1634- 671/1762)  16.7 E59-E36 N39-N51
      2 1119 1647 (N45W05 L=184)   281     41        1 (1096/1621-1141/1665)  56.3 W03-W07 N42-N46
      3  830 1632 (N42E19 L=160)   440     77        0 ( 795/1606- 865/1661)  38.7 E21-E16 N40-N44
      4 1230 1550 (N38W14 L=193)   277     51        1 (1208/1534-1261/1572) -36.9 W12-W16 N36-N39
      5  422 1365 (N21E45 L=134)   501    107        1 ( 397/1333- 452/1409) -45.0 E47-E41 N19-N24
      6  351 1280 (N15E49 L=130)  1262    236        4 ( 289/1232- 439/1365)   4.4 E55-E41 N12-N20
      7  538 1286 (N17E34 L=145)   235     50        2 ( 518/1259- 557/1313)  90.0 E35-E32 N15-N18
      8  339 1176 (N09E49 L=130)   859    259        1 ( 282/1140- 428/1240) -12.1 E54-E40 N06-N12
      9  196 1155 (N06E64 L=115)   476    113        0 ( 177/1112- 212/1206) -56.3 E67-E61 N03-N09
     10  937 1058 (N05E06 L=173)   201     88        0 ( 906/1021- 977/1091) -38.7 E07-E02 N02-N06
     11  608  940 (S04E26 L=153)   239     75        0 ( 591/ 899- 628/ 983)  78.7 E27-E24 S06-S01
     12  379  804 (S14E44 L=135)  1225    183        2 ( 350/ 727- 423/ 903)  79.7 E46-E40 S18-S07
     13 1858  843 (S07W65 L=244)   736    106        1 (1834/ 800-1876/ 881) -51.3 W61-W67 S09-S04
     14  960  781 (S12E03 L=176)   257     60        1 ( 934/ 744- 987/ 811)  63.4 E04-E01 S14-S10
     15 1418  581 (S24W29 L=208)   308     78        0 (1381/ 569-1459/ 597)  -9.5 W26-W32 S25-S23
     16 1224  431 (S36W18 L=197)   860    148        1 (1160/ 402-1299/ 463) -10.3 W12-W23 S38-S33
     17  829  405 (S40E13 L=166)  1391    327        2 ( 681/ 356- 975/ 439)  13.5 E26-E01 S44-S36
     18 1427  378 (S40W36 L=215)   335     58        0 (1408/ 355-1446/ 404)  71.6 W35-W37 S41-S37
  

    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