Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1702 (N46E48 L=128)  3589    511       11 ( 440/1621- 706/1757)  21.8 E56-E32 N38-N50
      2  864 1629 (N42E16 L=160)   451     93        1 ( 831/1602- 895/1668)  53.1 E18-E13 N40-N45
      3 1263 1553 (N38W16 L=192)   367     79        2 (1218/1524-1300/1584) -26.6 W13-W19 N36-N40
      4  444 1371 (N22E43 L=133)   776    192        3 ( 409/1328- 482/1429) -33.7 E47-E38 N19-N25
      5  257 1355 (N19E61 L=115)   617    126        1 ( 238/1294- 278/1415) -82.9 E62-E59 N15-N23
      6  403 1248 (N14E44 L=132)   983    184        5 ( 361/1202- 456/1288)  45.0 E47-E39 N11-N16
      7  380 1165 (N09E45 L=131)   749    249        2 ( 314/1135- 466/1207)  -8.7 E51-E37 N06-N10
      8 1883  853 (S07W69 L=245)   694     90        1 (1862/ 818-1901/ 887) -51.3 W65-W72 S09-S04
      9  419  782 (S15E41 L=135)  1025    203        3 ( 384/ 708- 445/ 855) -66.0 E43-E38 S19-S10
     10 1485  577 (S25W34 L=210)   464    188        0 (1421/ 550-1568/ 596) -14.0 W29-W41 S26-S23
     11 1262  428 (S36W21 L=196)   863    158        2 (1202/ 400-1342/ 464) -10.3 W15-W26 S38-S33
     12  910  410 (S39E06 L=169)   944    242        2 ( 802/ 372-1027/ 440)  19.4 E15-W02 S42-S35
  

    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