Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  715 1804 (N43E38 L= 27)   525     65        3 ( 689/1781- 745/1824)  56.3 E40-E35 N41-N45
      2  374 1474 (N16E51 L= 14)   310     66        0 ( 363/1434- 391/1505) -78.7 E52-E49 N13-N18
      3 1453 1290 (N15W24 L= 88)   308     50        0 (1421/1258-1479/1316) -26.6 W21-W25 N13-N16
      4  233 1091 (S08E56 L=  9)   832    195        5 ( 176/1044- 280/1129)  -4.8 E62-E50 S10-S05
      5  213 1026 (S12E57 L=  7)   474     96        2 ( 194/ 985- 236/1067)  78.7 E59-E54 S14-S09
      6  447 1022 (S11E37 L= 28)   269     69        0 ( 428/ 987- 463/1062)  76.0 E38-E35 S13-S08
      7 1642  849 (S08W42 L=107)  1232    283        3 (1551/ 812-1730/ 900)  12.1 W34-W48 S09-S04
      8  378  779 (S26E42 L= 22)   604    173        1 ( 348/ 718- 431/ 856)  76.0 E44-E38 S29-S21
      9 1070  738 (S22W06 L= 71)   188     64        0 (1045/ 702-1089/ 773)  78.7 W05-W07 S24-S19
     10 1471  588 (S26W36 L=101)   341     54        0 (1444/ 557-1494/ 626)  76.0 W34-W37 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