Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1079 1826 (N48E16 L=195)  1380    316        1 (1041/1739-1124/1906) -73.3 E19-E12 N39-N59
      2  380 1515 (N13E54 L=158)   454     85        0 ( 346/1491- 421/1534)  31.0 E56-E50 N11-N14
      3 1455 1407 (N24W20 L=232)   754    142        1 (1419/1348-1493/1458)  63.4 W18-W22 N20-N28
      4  364 1376 (N04E50 L=161)   565    107        0 ( 342/1316- 393/1431)  82.9 E51-E48 N00-N08
      5 1747 1290 (N24W44 L=255)  1185    140        2 (1703/1211-1782/1362) -53.1 W39-W47 N19-N28
      6  803  968 (S14E12 L=199)  1191    313        2 ( 690/ 832- 879/1117)  68.2 E17-E09 S23-S03
      7 1630 1046 (N06W37 L=248)   508    111        1 (1577/1008-1677/1100) -20.6 W32-W40 N04-N08
      8  329  928 (S23E45 L=166)   983    273        0 ( 238/ 841- 393/1016)  45.0 E53-E40 S28-S16
      9  434  774 (S31E36 L=176)  1489    392        1 ( 306/ 689- 568/ 892)  38.7 E47-E26 S37-S21
     10 1687  662 (S13W52 L=263)  3998    887        5 (1402/ 598-1891/ 729)  17.3 W30-W76 S19-S05
     11 1309  611 (S25W27 L=238)   327    107        1 (1257/ 588-1372/ 637)  23.3 W23-W29 S27-S22
     12  965  543 (S37W06 L=218)   239     73        0 ( 919/ 532-1006/ 554)  21.8 W03-W08 S38-S35
  

    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