Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1053 1593 (N33W02 L=158)   658    143        1 ( 972/1573-1134/1610)  -5.2 E03-W08 N31-N34
      2 1465 1159 (N05W27 L=183)   597     93        3 (1417/1105-1498/1208) -63.4 W24-W29 N01-N07
      3  273 1076 (N01E51 L=104)   738    141        3 ( 227/ 968- 310/1158) -65.6 E56-E47 S05-N06
      4  683 1099 (N01E21 L=135)   188     49        0 ( 653/1082- 713/1113) -26.6 E22-E18 N00-N02
      5  330  982 (S04E46 L=109)   582    163        1 ( 308/ 904- 351/1064)  77.5 E48-E44 S09-N00
      6  242 1015 (S02E55 L=101)   357     80        2 ( 218/ 970- 259/1058) -53.1 E57-E52 S04-N00
      7 1417  952 (S08W24 L=180)   586    197        1 (1354/ 850-1481/1021)  51.3 W20-W28 S13-S03
      8  197  981 (S04E60 L= 96)   401     91        2 ( 180/ 942- 222/1017)  59.0 E61-E56 S06-S01
      9  252  807 (S15E56 L=100)  3833    810       11 ( 159/ 710- 418/ 908)  -5.8 E65-E41 S20-S08
     10 1849  782 (S17W61 L=217)  3477    657        9 (1746/ 704-1939/ 869)  30.1 W51-W74 S21-S10
     11  416  807 (S15E41 L=115)   482     93        1 ( 387/ 772- 441/ 851)  51.3 E43-E38 S17-S12
     12 1680  754 (S19W45 L=201)   265     89        0 (1644/ 732-1705/ 794)  45.0 W42-W46 S20-S16
     13  313  667 (S23E52 L=104)   423     69        1 ( 289/ 627- 329/ 699)  76.0 E54-E50 S25-S21
     14  412  615 (S27E45 L=111)   550    156        0 ( 370/ 580- 468/ 642)  21.8 E49-E39 S29-S25
     15  292  584 (S29E57 L= 99)   461    142        3 ( 252/ 564- 321/ 612) -18.4 E60-E53 S29-S26
     16  599  522 (S34E32 L=124)   773    216        1 ( 579/ 438- 633/ 611)  63.4 E35-E27 S40-S28
     17 1583  557 (S32W41 L=197)   365     96        1 (1543/ 536-1620/ 592)  38.7 W38-W43 S33-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