Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  804 1632 (N30E31 L= 98)   352     84        1 ( 761/1599- 833/1649)  21.8 E34-E28 N27-N31
      2  545 1523 (N17E44 L= 85)   654    163        2 ( 483/1486- 591/1555)  12.5 E48-E39 N14-N18
      3  804 1441 (N17E23 L=106)   628    186        2 ( 735/1394- 882/1463)  32.0 E27-E18 N13-N20
      4 1402 1379 (N28W16 L=145)   543     96        1 (1370/1331-1423/1427) -71.6 W13-W17 N24-N30
      5  750 1338 (N10E23 L=106)   267     86        0 ( 703/1327- 789/1351)  21.8 E26-E21 N08-N11
      6  575 1147 (S05E30 L= 99)   363     84        0 ( 550/1111- 606/1181)  68.2 E30-E28 S07-S02
      7  552  859 (S22E26 L=103)  1162    355        4 ( 442/ 761- 648/ 999)  58.6 E34-E20 S30-S12
      8  277  756 (S35E52 L= 77)   381     31        0 ( 263/ 729- 296/ 780)  71.6 E54-E49 S36-S33
      9 1379  520 (S24W36 L=166)   534    105        1 (1343/ 492-1416/ 570)  90.0 W34-W40 S25-S20
     10  720  442 (S46E07 L=122)  1385    249        4 ( 607/ 420- 842/ 468)  19.3 E17-W03 S50-S41
     11 1120  203 (S54W43 L=172)  1302    207        1 (1069/ 136-1177/ 287) -46.6 W36-W54 S62-S44
  

    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