Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1087 1745 (N40E22 L= 61)   320     55        0 (1063/1728-1118/1764) -11.3 E23-E18 N38-N41
      2  507 1743 (N23E65 L= 18)   789     79        0 ( 458/1719- 536/1761)   0.0 E73-E59 N21-N24
      3  804 1682 (N26E35 L= 47)   274     56        0 ( 775/1666- 841/1694)  26.6 E37-E32 N24-N27
      4 1312 1497 (N30W03 L= 86)  1019    213        1 (1232/1406-1366/1574) -28.6 E03-W08 N24-N33
      5  983 1411 (N14E13 L= 69)   326    142        0 ( 964/1348-1007/1474) -76.0 E15-E11 N10-N18
      6 1106 1316 (N12E04 L= 79)   310    115        0 (1072/1264-1130/1355) -51.3 E06-E01 N09-N14
      7 1653 1205 (N23W34 L=116)   583    106        1 (1612/1151-1678/1240) -31.0 W29-W36 N20-N25
      8 1342 1177 (N11W13 L= 96)   686    244        1 (1232/1095-1386/1237) -14.0 W05-W17 N07-N13
      9 1627 1078 (N15W34 L=117)   370     68        1 (1603/1043-1646/1117)  90.0 W32-W35 N12-N17
     10 1799  991 (N16W51 L=134)  1233    233        2 (1709/ 925-1867/1042)   3.0 W41-W60 N13-N17
     11 1720  924 (N09W46 L=128)   480    124        0 (1698/ 876-1754/ 987) -35.5 W42-W50 N07-N12
     12 1233  501 (S28W30 L=113)   445    110        0 (1179/ 485-1284/ 520)  33.7 W27-W33 S30-S26
     13  644  468 (S51E09 L= 74)   398     70        0 ( 622/ 446- 675/ 501) -26.6 E11-E05 S52-S48
     14  779  482 (S45W02 L= 85)   236     40        0 ( 762/ 463- 803/ 501)  76.0 W00-W03 S46-S42
     15 1596  387 (S22W65 L=148)   461    112        2 (1559/ 345-1622/ 414) -18.4 W61-W72 S23-S20
     16 1216  353 (S38W39 L=122)   540    155        0 (1157/ 306-1275/ 394)  77.5 W37-W40 S42-S33
     17 1372  343 (S32W51 L=133)   693    194        0 (1298/ 297-1457/ 377)  63.4 W47-W54 S37-S27
     18 1456  322 (S31W59 L=142)   399     68        0 (1430/ 292-1482/ 353) -59.0 W57-W61 S32-S27
  

    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