Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1421 1376 (N13W33 L=241)   265     80        0 (1387/1350-1449/1404)  11.3 W29-W35 N12-N14
      2 1626 1276 (N02W43 L=251)   442    140        0 (1582/1216-1671/1331)  16.7 W37-W48 N00-N03
      3 1177 1158 (N08W12 L=220)   919    236        3 (1106/1088-1245/1211)   0.0 W06-W16 N04-N11
      4  138 1141 (N32E67 L=141)   641     73        0 ( 109/1110- 163/1168) -20.0 E74-E63 N29-N33
      5  565 1085 (N19E26 L=183)   248    117        1 ( 497/1076- 613/1106) -14.0 E30-E22 N17-N20
      6  389 1008 (N19E40 L=168)   370     93        1 ( 367/ 957- 409/1049)  38.7 E42-E37 N16-N20
      7  785  854 (N01E17 L=191)   365    128        1 ( 744/ 792- 821/ 902) -80.5 E18-E16 S02-N04
      8 1467  741 (S23W19 L=227)  1277    313        4 (1394/ 613-1527/ 870)  38.2 W11-W25 S28-S16
      9 1069  817 (S08E03 L=205)   230     72        1 (1036/ 791-1092/ 849)  26.6 E05-E01 S08-S06
     10  696  711 (S05E27 L=182)   422    195        0 ( 664/ 633- 765/ 755)  26.6 E29-E21 S08-S02
     11 1875  728 (S38W59 L=267)   977     78        3 (1838/ 700-1908/ 755) -18.4 W52-W67 S40-S35
     12  972  670 (S14E12 L=196)   374    161        1 ( 928/ 614-1032/ 702)  26.6 E14-E08 S16-S11
     13  307  555 (S05E62 L=146)  4814    378        8 ( 217/ 381- 375/ 685)  53.6 E83-E55 S16-N03
     14  814  559 (S16E25 L=184)   247     20        0 ( 790/ 540- 836/ 580)  14.0 E26-E22 S17-S14
     15 1266  436 (S37E02 L=206)   284     82        0 (1239/ 404-1297/ 469)  26.6 E04-W00 S38-S34
     16  562  333 (S24E54 L=154)  1348    174        3 ( 492/ 294- 638/ 386)   8.5 E64-E44 S25-S22
  

    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