Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1479 (N17E33 L= 85)   658    155        3 ( 613/1437- 737/1520)   0.0 E37-E27 N14-N19
      2  744 1428 (N15E27 L= 91)   279     74        2 ( 719/1398- 775/1460)  63.4 E28-E25 N13-N17
      3  615 1419 (N12E35 L= 83)   336     76        1 ( 592/1384- 643/1449)  90.0 E36-E33 N09-N14
      4  979 1392 (N18E11 L=107)   299    111        1 ( 929/1371-1030/1414)  15.9 E14-E07 N16-N20
      5  725  811 (S21E14 L=104)   343    101        0 ( 698/ 770- 763/ 839)  90.0 E15-E11 S23-S18
      6  628  729 (S28E19 L= 99)   600    148        1 ( 564/ 699- 694/ 758)  33.7 E23-E14 S31-S25
      7  393  700 (S36E39 L= 79)   406     63        1 ( 373/ 663- 408/ 728) -53.1 E41-E37 S38-S33
      8 1500  457 (S24W49 L=167)   645     56        2 (1471/ 425-1539/ 487)   0.0 W46-W53 S26-S22
      9  573  450 (S50E22 L= 96)   417     71        0 ( 536/ 435- 609/ 465)  18.4 E26-E17 S51-S48
     10  810  390 (S48W02 L=120)   826    152        1 ( 732/ 368- 882/ 414)  18.4 E05-W08 S50-S45
     11 1378  346 (S35W48 L=165)   385     87        0 (1348/ 317-1418/ 371)  90.0 W46-W49 S37-S31
  

    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