Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1047 1575 (N40W03 L=223)   826    220        1 ( 965/1515-1137/1640) -36.0 E03-W09 N35-N45
      2 1612 1463 (N29W47 L=267)   819    158        2 (1538/1425-1683/1498) -17.1 W41-W54 N26-N32
      3  196 1262 (N18E67 L=152)  2433    158        8 ( 154/1189- 232/1354) -78.7 E73-E62 N13-N23
      4 1588 1263 (N17W39 L=259)   318     90        0 (1542/1247-1628/1277)   0.0 W35-W42 N15-N17
      5  679 1202 (N15E22 L=197)   304     79        2 ( 663/1163- 698/1248) -80.5 E23-E20 N12-N18
      6  272 1151 (N11E55 L=165)   651     95        2 ( 237/1103- 313/1187)  39.8 E58-E50 N08-N13
      7 1637 1124 (N08W41 L=261)   384    130        1 (1615/1060-1656/1181)  60.3 W39-W43 N04-N11
      8  478  854 (S06E36 L=183)  1149    264        3 ( 392/ 767- 574/ 951) -37.6 E42-E29 S11-N00
      9  283  836 (S08E54 L=166)  1044    158        4 ( 261/ 761- 302/ 902) -83.7 E55-E51 S12-S03
     10  334  790 (S11E49 L=170)   379     59        2 ( 311/ 755- 349/ 825)  63.4 E51-E47 S12-S08
     11  957  748 (S13E05 L=215)   425    127        1 ( 915/ 704- 981/ 805)  56.3 E07-E03 S15-S09
     12  323  660 (S19E54 L=166)  4201    461       11 ( 180/ 531- 457/ 794) -33.7 E68-E43 S27-S11
     13  553  608 (S22E34 L=186)  1178    211        1 ( 504/ 509- 610/ 689)  49.8 E39-E28 S29-S16
     14 1103  578 (S25W04 L=224)   747    175        0 (1059/ 499-1154/ 642)  59.0 W01-W07 S30-S20
     15 1209  575 (S25W12 L=231)   294    120        0 (1157/ 536-1260/ 594)  26.6 W08-W15 S27-S23
     16  508  378 (S40E48 L=171)   471     64        0 ( 497/ 339- 519/ 417)  45.0 E52-E44 S43-S36
     17  912  352 (S42E11 L=209)  1597    259        3 ( 803/ 315-1014/ 405)  15.5 E20-E02 S45-S37
     18 1265  177 (S63W30 L=250)   876     74        1 (1247/ 136-1283/ 212) -34.4 W24-W43 S70-S57
  

    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