Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  725 1839 (N58E52 L=309)  1764    104        2 ( 681/1801- 777/1870)  21.8 E62-E44 N53-N61
      2  857 1612 (N39E19 L=342)   992    208        3 ( 778/1569- 957/1652) -13.2 E26-E09 N36-N41
      3 1099 1448 (N30W02 L=  2)   293     65        1 (1079/1428-1130/1480) -36.9 W00-W04 N28-N31
      4  235 1417 (N20E71 L=289)  1445    104        4 ( 206/1368- 264/1462) -14.7 E86-E64 N16-N22
      5 1489 1401 (N29W31 L= 31)   260     78        0 (1465/1370-1517/1434) -53.1 W29-W33 N27-N31
      6  626 1338 (N18E29 L=331)   354    135        0 ( 574/1286- 660/1376) -32.0 E33-E25 N15-N20
      7 1739 1292 (N23W53 L= 53)   489    123        2 (1701/1252-1776/1335)  29.1 W47-W58 N20-N26
      8 1711 1181 (N16W48 L= 48)   348     46        1 (1695/1146-1733/1207) -90.0 W46-W50 N13-N17
      9 1141 1084 (N06W07 L=  7)   353     60        1 (1106/1052-1174/1104)  14.0 W04-W08 N04-N07
     10 1705  647 (S17W54 L= 54)   454     44        1 (1681/ 627-1731/ 672)  71.6 W51-W57 S18-S15
     11 1464  558 (S25W35 L= 36)   424    122        0 (1411/ 541-1518/ 586) -11.3 W30-W40 S25-S23
  

    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