Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1053 1622 (N39W22 L=161)   879    207        1 (1012/1539-1107/1690)  33.7 W16-W29 N33-N41
      2 1468 1594 (N24W47 L=185)  1045    157        2 (1403/1550-1539/1620) -45.0 W41-W50 N19-N26
      3  513 1241 (N31E26 L=112)  1030    351        0 ( 375/1210- 675/1277) -22.2 E37-E15 N26-N36
      4 1522 1155 (S01W32 L=171)   978    231        2 (1417/1119-1612/1187) -17.1 W24-W37 S05-N02
      5 1470 1074 (S04W26 L=165)   310    117        0 (1434/1027-1522/1099)   8.1 W23-W30 S05-S02
      6 1683  977 (S16W38 L=177)   741    107        3 (1636/ 932-1726/1016)  23.2 W34-W42 S17-S14
      7 1414  827 (S16W17 L=155)   707    210        2 (1390/ 748-1456/ 947)  81.3 W14-W19 S21-S08
      8 1083  683 (S15E06 L=133)   711    164        1 ( 990/ 644-1161/ 713) -41.2 E10-E02 S19-S11
      9 1178  577 (S24E04 L=135)   367     80        0 (1162/ 534-1193/ 621)  71.6 E04-E01 S27-S21
     10  754  494 (S18E31 L=107)   264    123        1 ( 696/ 473- 806/ 513)  -7.1 E35-E27 S18-S16
     11  658  439 (S19E40 L= 99)   330     79        0 ( 615/ 410- 680/ 467) -63.4 E41-E38 S20-S16
     12  777  399 (S24E34 L=104)   521    111        3 ( 747/ 362- 808/ 439) -90.0 E36-E32 S27-S21
     13 1030  214 (S46E30 L=108)  1048    151        2 (1000/ 172-1060/ 256)  28.6 E36-E25 S49-S42
  

    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