Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  533 1422 (N03E40 L=131)   332     61        2 ( 515/1393- 559/1447) -56.3 E42-E38 N01-N04
      2  854 1359 (N07E18 L=153)   241     84        0 ( 810/1335- 900/1378)  51.3 E20-E16 N04-N09
      3 1337 1327 (N19W09 L=181)   371    147        1 (1269/1283-1415/1362)   0.0 W05-W15 N16-N19
      4  293 1232 (S12E53 L=118)   327    158        0 ( 252/1194- 331/1276) -16.7 E58-E48 S13-S10
      5 1263 1112 (N04W11 L=182)   316    104        2 (1212/1095-1319/1130)  15.9 W07-W14 N02-N05
      6  463  880 (S29E32 L=140)   973    260        1 ( 356/ 823- 551/ 936)  42.7 E39-E26 S35-S23
      7  889  856 (S20E03 L=168)   190     28        1 ( 865/ 836- 907/ 873)   0.0 E04-E01 S20-S18
      8 1647  698 (S06W48 L=219)   630    199        2 (1581/ 661-1708/ 720)   4.8 W42-W54 S07-S04
      9 1680  562 (S11W58 L=229)   472    162        2 (1634/ 538-1739/ 578)  29.1 W53-W63 S13-S08
     10 1022  502 (S36W17 L=188)  1098    251        4 ( 919/ 437-1113/ 575)  -6.7 W08-W25 S38-S33
     11 1342  545 (S24W34 L=205)   323    101        0 (1280/ 523-1388/ 569)  50.2 W31-W36 S26-S20
     12 1290  449 (S30W36 L=207)   332    103        1 (1258/ 417-1326/ 491) -14.0 W31-W39 S31-S29
  

    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