Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  987 1892 (N61W48 L= 33)   456     81        0 ( 953/1858-1021/1918)   3.2 W38-W56 N60-N62
      2  660 1709 (N60E05 L=341)   371     61        0 ( 623/1689- 691/1723)  -6.3 E09-W00 N58-N62
      3 1721 1513 (N12W63 L= 49)  6873      0       14 (1533/1343-1827/1633) -80.1 W53-W78 N00-N23
      4 1182 1554 (N31W27 L= 12)   476    146        0 (1117/1523-1245/1593)   0.0 W21-W32 N30-N32
      5  677 1454 (N40E10 L=335)   810    200        0 ( 591/1423- 758/1477) -13.0 E17-E04 N36-N43
      6 1265 1453 (N23W27 L= 12)   334    108        0 (1217/1436-1323/1473)  -7.1 W23-W31 N21-N24
      7 1385 1443 (N19W34 L= 19)   238     94        0 (1349/1413-1419/1462) -11.3 W30-W35 N17-N20
      8  305 1362 (N44E48 L=298)  1001    161        1 ( 272/1306- 351/1435) -55.0 E53-E41 N39-N49
      9 1257 1307 (N15W21 L=  6)   251    104        0 (1243/1253-1272/1355)  59.0 W19-W22 N12-N17
     10 1846 1281 (S04W63 L= 49)  1322    110        0 (1787/1247-1887/1328)   0.0 W55-W72 S05-S02
     11  852 1232 (N22E04 L=341)   595    199        0 ( 787/1182- 924/1263) -29.7 E09-E01 N17-N24
     12 1691 1166 (S05W44 L= 29)   620    154        0 (1667/1093-1725/1230)  40.6 W40-W48 S08-S02
     13  894 1114 (N15E05 L=340)   352     98        0 ( 867/1060- 922/1166)  38.7 E07-E02 N11-N17
     14  213 1080 (N28E57 L=288)   402     96        0 ( 201/1025- 226/1131) -81.9 E58-E55 N24-N31
     15 1147  622 (S19E04 L=342)   634    207        0 (1090/ 527-1185/ 711)  45.0 E09-W00 S23-S14
     16  597  583 (S08E37 L=308)   933    210        0 ( 570/ 496- 627/ 699)  47.7 E42-E32 S12-S01
     17 1710  543 (S45W37 L= 23)  1512    200        2 (1621/ 488-1796/ 601)  -2.6 W27-W50 S46-S42
     18  772  249 (S34E45 L=301)  2173    321        2 ( 641/ 166- 885/ 332)   8.7 E69-E30 S38-S30
  

    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