Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1299 1393 (N24W07 L= 36)   836    252        3 (1238/1280-1370/1493) -37.6 W00-W13 N19-N29
      2  842 1307 (N07E18 L= 12)   380    150        0 ( 778/1279- 925/1338)  45.0 E20-E13 N04-N11
      3  488  900 (S24E30 L=359)   810    195        2 ( 410/ 862- 568/ 953)  36.0 E36-E25 S27-S19
      4 1125  841 (S12W10 L= 40)   233    120        0 (1074/ 819-1193/ 864)  50.2 W08-W13 S14-S08
      5 1161  794 (S13W14 L= 44)   250     81        0 (1124/ 782-1207/ 811)  31.0 W11-W16 S14-S11
      6  876  724 (S25E01 L= 29)   314     91        1 ( 826/ 702- 931/ 755)  18.4 E04-W02 S25-S23
      7  457  720 (S36E31 L=359)   259     66        1 ( 431/ 698- 480/ 743)  36.9 E32-E28 S38-S34
      8  543  632 (S40E22 L=  8)  1233    195        2 ( 466/ 551- 624/ 695) -14.9 E29-E14 S43-S36
      9 1452  614 (S15W38 L= 68)   521    142        0 (1419/ 549-1495/ 680) -69.4 W36-W40 S18-S10
     10  770  567 (S38E03 L= 27)   392     95        0 ( 723/ 551- 825/ 582)  18.4 E06-W00 S39-S36
     11 1301  547 (S23W31 L= 61)   596    245        0 (1196/ 531-1421/ 560)  21.8 W24-W39 S26-S20
     12 1596  339 (S26W65 L= 95)   577     97        1 (1563/ 320-1627/ 358)  45.0 W62-W67 S27-S23
  

    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