Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  385 1083 (N23E42 L= 62)   268     50        0 ( 369/1054- 400/1111)  56.3 E43-E40 N21-N25
      2 1700  987 (S12W42 L=146)   329     38        0 (1670/ 972-1728/1002) -21.8 W39-W44 S13-S11
      3  413  845 (N09E43 L= 61)   593    143        2 ( 386/ 803- 444/ 896)  68.2 E45-E40 N06-N11
      4 1187  620 (S20W00 L=104)  1516    503        1 (1017/ 464-1330/ 782)  26.6 E14-W12 S26-S13
      5 1428  617 (S27W16 L=119)   408    172        0 (1376/ 564-1454/ 698)  45.0 W11-W18 S29-S22
      6  882  567 (S17E20 L= 84)   226     91        0 ( 858/ 534- 917/ 607)  38.7 E21-E16 S18-S14
      7  896  420 (S27E24 L= 80)  1387    312        2 ( 832/ 373-1014/ 529) -77.0 E27-E16 S31-S18
      8 1241  427 (S35E01 L=102)   251     55        1 (1214/ 403-1269/ 452)  14.0 E03-W01 S36-S34
      9  587  261 (S32E60 L= 43)   689     84        1 ( 558/ 221- 607/ 292)  45.0 E66-E55 S36-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